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Abstract 


A  tunable  diode  laser  is  used  to  investigate  the  v  ^  vibrational  band  of  Sulfuryl  Fluoride  in  the 

11  g  region.  Although  the  absorption  spectrum  is  doppler-limited,  the  resolution  is  high 

enough  to  observe  previously  unresolved  lines  and  therefor  determine  the  band  center  and  other 
excit^  state  molecular  constants. 

The  band  constants  are  determined  by  a  least  squares  fit  computer  program  using  an  S 
Reduced  hamiltonian  to  model  the  molecular  vibration-rotation.  The  molecule  is  slightly 
aspherical  and  requires  curve  fitting  to  thousands  of  lines. 

The  S  Reduced  hamiltonian,  as  well  as  other  pertinent  theoretical  topics  are  discussed. 
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Chapter  1 


Introduction 

The  advent  of  infrared  lasers  h?s  led  to  unprecedented  high-resolution  spectroscopy  in 
the  middle  IR  region.  The  importance  of  the  laser  in  spectroscopy  is  evident  in  the  study  of  the 
molecule  Sulfuryl  Fluoride  (SF202)-  The  availability  of  high-rcsolufion  data  allows  us  to 
observe  previously  unresolved  structure,  and  for  the  first  time,  observe  such  effects  as 
asymmetric  K  splitting  and  the  impact  of  nuclear  spin  statistics  in  this  molecule,  as  well  as 

determine  band  centers  and  upper  state  molecular  constants  to  a  high  degree  of  accuracy.  In  this 
experiment  a  tunable  diode  laser  (TDD  was  used  to  map  the  spectrum  of  SF2O2  in  the  11  |i 

range. 


1.1  Previous  Studies  of  SF2O2 


Microwave  studies  of  SF2C2  have  revealed  seven  lines  in  the  20  to  30  thousand  Mhz 
region.  Six  of  these  lines  were  assigned  to  the  rotational  transitions  J  =  1  — >  2  and  J  =  2— >3 
of  ^^F202  and  one  line  to  the  rotational  transition  J  =  1  — >  2  of  the  ^^SF202  isotope.  With 
this  data,  and  assuming  C2v  symmetry,  the  structure  of  SF2O2  was  determined  and  the  ground 
state  molecular  constants  estimated.  .Microwave  studies  have  also  determined  that  SF2O2  is 
slightly  asymmetrical.  However,  due  to  the  smallness  of  deviation  from  sphericity  of  SF2O2/ 
the  b  and  c  axes  (principle  axes  of  intermediate  and  greatest  moments  of  inertia  respectively) 
could  not  be  unambiguously  assigned  to  their  respective  molecular  axes.  In  addition  the 
microwave  data  obtained  was  not  detailed  enough  to  determine  molecular  distortion  parameters 
(14). 
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There  has  also  been  infrared  spectroscopy  performed  on  SF2O2  in  which  band  centers  of 
fundamental  and  combination  vibrations  have  been  recorded,  but  the  transition  lines  themselves 


were  unresolved  [11]  [12]  (see  Figure  1). 


Wavenumbers  in  CM 


Figure  1~  Infrared  Spectrum  of  SF2O2  between  250  cm'^  and  4000  cm*V  From  (12]  (p.  571]. 

This  spectrum,  recorded  in  1960,  fails  to  resolve  any  but  the  coarsest  of  features.  This 
means  that,  in  the  Vg  vibrational  band  for  example  (depicted  between  the  arrows  in  Figure  1), 

there  are  approximately  250  transitions  and  thousands  of  strong  lines  that  he  entirely  below  the 
instrumental  resolution  of  the  times.  The  power  of  our  technique  lies  in  our  ability  to  take  a 
segment  of  this  previously  recorded  spectrum  and  resolve  completely  all  the  transitions  and 
strong  lines.  In  particular,  our  work  resolves  transitions  over  a  wide  range  extending  from  J  =  30  of 
the  R-branch  to  J  =  16  of  the  P-branch  of  the  Vg  vibrational  band.  The  significance  and 


7 


19.63 


869.65 


wavenumbsr(cm-l) 

Rgure  2  -  Experimental  Spectra.  The  Resolution  is  about  60  MHz 

J=7-6 


v/avenumb8r{cm-1 ) 


resolution  of  our  IR  data  (see  Figure  2  and  results  in  Chapter  5)  is  apparent  when  compared  to 
the  previously  recorded  IR  spectra  shown  in  figure  1. 

In  addition  to  our  achieved  resolution  (w’hich  is  doppler-limited  to  about  60  Mhz  at 
room  temperature),  we  have  identified  thousands  of  lines  using  a  computer  generated  least 
squares  fit.  Using  the  fit,  we  have  been  able  to  clarify  the  assignments  of  the  b  and  c  axes  to  the 
appropriate  SF2O2  molecular  axes  and  determine  the  excited  state  inertial  constants  A,  B  and  C 
the  Vg  band  center,  and  the  excited  state  distortion  parameters.  Xone  of  this  information  had 
been  previously  available  except  for  the  band  center,  which  had  been  estimated  (to  the 

nearest  wavenumber)  at  885  cm"^.  To  illustrate,  some  of  the  measured  V'alues  from  this 
experiment  are  given  in  Table  1. 


1.2  Thesis  Organization 

The  purpose  of  this  thesis  is  to  identify  completely  the  P  and  R-branch  lines  covering  a 
15  cm’^  range  around  the  Q-branch  of  the  V  g  vibrational  band  of  SF2O2,  and  determine  the 

upper  state  rotational  and  centrifugal  band  constants  as  well  as  the  Vg  vibrational  band  center. 

Normally  one  might  expect  to  perform  an  experiment  to  check  the  validity  of 
theoretically  derived  parameters.  In  this  experiment  we  have  the  opposite  situation.  Here  we 
are  empirically  deriving  values  from  observed  spectral  frequencies  which  arc  then  placed  into 
the  theoretical  equation.  This  technique  is  described  as  a  reverse  stochastic  process  in  which 
asymmetric  rotor  quantum  numbers  arc  assigned  to  observ'Cd  transitions  using  av'ailablc  theory. 
Once  quantum  numbers  are  assigned,  band  constants  are  determined  by  curve  fitting  using  a 
computerized  least  squares  routinelSOj.  The  curve  to  be  fitted  is  the  S-reduced  hamiitonian 
derived  in  .Appendix  A. 

There  are  several  steps  in  identifying  and  fitting  the  SF2O2  spectra.  These  steps  are 
covered  in  the  various  chapters  of  this  thesis. 

Chapter  2  outlines  the  basic  theory  of  asymmetric  molecules  and  establishes  the 
notations  and  conventions  used  in  the  thesis. 

Chapter  3  explains  the  experimental  'ctup  and  the  procedure  used  to  collect  the  spectra. 
.Also  covered  is  the  application  of  the  basic  theory  to  the  interpretation  and  analysis  of  the 
experimental  spectra. 

Chapter  4  covers  error  analysis. 

The  final  results  of  the  experiment  are  given  in  Chapter  5.  Over  46vX)  lines  were  fitted 
by  the  least  squares  fitting  routinefSO].  Instead  of  listing  all  46(K)  line5,sevcral  representative  P 
and  R-branch  transitions  are  presented  in  graphical  form. 

The  development  of  the  S  Re  duoKi  hamiitonian  and  the  wave  equation  solution  for  the 
symmetric  rotor  are  the  subjects  of  Appendix  .A  and  B  respectively 
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Chapter  2 
Basic  Theory 

The  original  contribution  of  this  thesis  begins  in  Chapter  3  on  page  36.  But  first,  the 
basic  theory  of  asymmetric  vib-rotors  is  given  in  this  chapter,  as  follows. 

2.1  The  Rigid  Rotor-Harmonic  Oscillator  Approximation^ 


To  begin  determining  the  band  constants  and  fitting  the  observ-ed  spcchra,  it  is  necessary 
to  use  any  good  approximation  that  is  available.  Oftentimes  approximations  can  lead  to 
reliable  estimates  of  such  quantities  as  band  constants  or  transition  spadngs,  whidt  can  then  be 
inserted  into  a  least-squares  program  to  begin  the  fitting  procedure  from  which  accurate 
molecular  parameters  are  determined.  Estimated  quantities  can  also  ser\'e  as  a  reliable  means  of 
calculating  errors  and  ensuring  the  fitting  process  stays  accurate. 

The  best  approximation  is  called  the  rigid  rotor-harmonic  oscillator  approximation. 
This  approximation  begins  with  the  hamiltonian  below'  (see  .Appendix  A) 


H  = 


Cj-Px)’  (‘’y-Pv)’ 


Cz-Pj  ■  I 


2! 


XX 


2Jvv 


2! 


zz 


-iZpl-v 

k 


(1) 


Px,  Py  and  Pz  components  of  the  total  angular  momentum,  pxx  Py  and  pz  are  components  of  the 

internal  angular  momentum,  pk  is  the  momentum  associated  with  the  kth  normal  vibrational 
mode,  and  V  represents  the  overall  molecular  potential.  The  first  three  terms  represent  the 


^The  theory  of  Section  2.1  is  taken  from  the  following  References.{31  (pp.  7-91;  12)  Ipp.  ©-541;  (251 
fpp.  266-267);  [7)  [pp.  2S-34I. 
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rotational  energy  and  the  interaction  of  the  rotational  and  internal  an^lar  nxinK^tuni.  The  last 
two  terms  are  pure  functions  of  the  molecular  vibrations. 


One  can  make  further  approximations  bj  neglecting  the  interna!  (vibrational)  angular 
momentum  and  approximating  the  vibrational  part  cf  the  hamiltonian  as  a  hariiKmic  osdilator 
potential  which  can  be  solved  exactly.  In  this  wav  the  rotational  and  vibrational  emsgies  are 
completely  separable  and  arc  represented  by  the  energy  levels  of  a  rigid  rotor  and  a  harmonic 

oscillator,  fn  this  experiment  the  energy  solution  to  the  vibrational  part  of  the  hamiltonian  is 
represented  by  the  band  center  .  The  remaining  (rotatiortal)  part  of  the  hamiltonian  b  the 

problem  we  raw  solve. 

The  energy  levels  of  a  rigid  rotor  arc  found  by  solving  the  SchrcKxiin^jr  equation: 


where,  with  all  the  above  approximations,  wo  nov.-  have 


(2) 


(3) 


where  J  now  symbolizes  the  total  angular  momentum  and  I  designates  tra  ngid  moment  of  irsstia 
about  the  indicated  axis.  An  asymmetric  lop  is  dtrfirad  by  the  relation 

i  ^  yy  ^  ^  7z  (4) 

Thb  relation  removes  tra  K-degcncracy  of  the  symoKitrsc  rotor  and.  as  a  result,  complKatcs  the 
rigid  rotor  spectrum.  Sirac  the  total  angubr  momentum  |  and  lU  projection  .M  on  a  spare-fiMSd 
axis  are  omstants  of  the  iralion.  iNiy  arc  ^x?d  qu,>nlum  nurraers  ar^J  can  be  used  to  speafv  iN? 
state.  However,  the  projection  of  J  on  a  molecule  lixcd  avis  is  no  loiter  constant  and  Iteefore  we 
can  not  use  K  to  speafy  il«:  state.  Instead  we  create  a  parameter  tt«  indicate  vznmis  cearces  el 
deviation  from  the  symmetric  prolate  ta.vis  of  symirairy  is  tN:  axis  of  least  mon^t  of  irartia) 


and  symmetric  oblate  (axis  of  symmetry  is  the  axis  of  greatest  moment  of  inertia)  limits.  The 
most  common  parameter  is  the  asymmetry  parameter 


2B-A- 

=  “-C'  (5) 


which  varies  from  -1  for  a  prolate  symmetric  top  (B=C)  to  1  for  an  oblate  symmetric  top  (B=A). 
The  quantities  A,  B  and  C  are  constants  inversely  proportional  to  the  least,  intermediate  and 
greatest  moments  of  inertia  respectively.  They  are  defined  in  Appendix  B  (see  p.  117).  Energy 

levels  are  specified  by  the  symbol  J  k  ^  •  This  level  represents  a  rotational  level  of  J  which 

in  the  limiting  case  connects  a  prolate  symmetric  top  level  of  K  with  an  oblate  symmetric  top 
level  of  (see  Figure  3). 

We  now  expand  the  asymmetric  top  eigenfunctions  in  combinations  of  symmetric  top 
basis  functions 


V,M(e,o,x)=2;A,lJKM> 

k 

Where  [JKM)  represents  the  symmetric  top  wave  function  of  Equation  (B-15)  (see  p.  117). 
Substituting  Equation  (6)  into  Equation  (2)  yields 

Xa^hjjkm)  =  e£aJjkm> 

k  k 
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The  energies  are  the  roots  of  the  (2J  +  1)  X  (2J  +  1)  secular  determinant 

-  E5,,j,.|  =  0  (8) 

where 

H^^.=  aK'M|H|,|JKM>  (9) 

At  this  point  it  is  helpful  to  investigate  the  symmetry  properties  of  the  wave  functions. 
For  this  we  resort  to  group  theory.  In  group  theoretical  nomenclature,  all  rigid  rotor  wave 
functions  belong  to  representations  called  the  external  rotation  group.  These  representations  are 
characterized  by  J  and  and  correspond  to  an  infinitely  large  amount  of  possible  rotations  of  the 
space-fixed  system  that  leave  the  energy  invariant.  The  symmetric  rotor,  described  by  the 
body-fixed  axes,  belongs  to  the  internal  rotation  group  D  o^/  characterized  by  the  quantum 
numbers  J  and  K.  For  the  asymmetric  rotor,  the  hamiltonian  has  the  character  that  a  rotation  of 
180°  about  any  axis  leaves  the  hamiltonian  unchanged  (see  Equation  (3)].  The  asymmetric  top 
is  therefore  classified  as  a  member  of  the  point  group  D  ,  ,  also  known  as  the  Viergruppe  tV) 

group.  This  group  is  characterized  by  twofold  rotations  about  the  a,  b  and  c  axes  denoted  by 
^2(a)'  *"2(b)^”'^^2(c)  and  the  identity  E  (see  Table  2). 


Table  2 

Group  Characters  for  the  D2(V)  Group.  From  [25]  (p.  265]. 


Group  Operations 

(Quantum 

Numbers 

K.,K, 

E 

C2(c) 

C2(b) 

1 

1 

1 

1 

ee 

1 

1 

-1 

-1 

oe 

1 

-1 

1 

-1 

00 

1 

.1 

-1 

1 

eo 
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Here  a,  b  and  c  represent  principle  molecular  axes  of  least,  intermediate  and  greatest 


moments  of  inertia  respectively.  They  correspond  to  the  inertial  constants  A,  B  and  C  defined  in 
Equation  (B-l(  '  (see  p.  117).  The  letters  "e"  and  "o"  in  Table  2  stand  for  even  and  odd 
respectively. 

Now  we  examine  how  the  |JKM)  transform  under  the  operations  of  the  V  group.  First 
we  note  that  for  SF2O2,  K  =  --5  [14],  which  approaches  the  symmetric  prolate  limit.  Therefore 

it  is  natural  to  identify  the  X,  Y,  Z  axes  with  the  b,c,a  axes  respectively  (see  Table  3). 


Table  3 

Possible  Identification  of  Principle  Axes  a,  b,  c  with  Body-fixed  Coordinates  x,  y,  z. 

From  [25]  |p.  268). 


Body-fixed 

coordinates 

Representation 

la 

me 

X 

b 

c 

a 

y 

c 

a 

b 

z 

a 

b 

c 

^most  natural  limit  of  a  prolate  symmetric  rotor 

^  most  natural  for  a  rotor  intermediate  between  the  two  limiting  symmetric  rotor  cases 
^  most  natural  limit  of  an  oblate  symmetric  rotor 


The  effect  of  the  V  group  operations  on  the  |JKM)  is  determined  by  their  effect  on  the  Euler 
angles.  We  use  the  direction  cosine  matrix  ^  p(0  ^  of  Table  B-1  to  determine  the  effects  of 


each  of  the  four  operations  on  the  Euler  angles.  The  results  are  summarized  in  Table  4. 


Table  4 

Cartisian  Coordinate  and  Euler  Angle  Transformations  Under  Operations  of  the  D2(V)  group. 

From  [25]  (p.  269], 


Group  Operation 

Cartisian  Coordinate 
Transformation 

Euler  Angle 
Transformation 

E 

X  -^x 

y-^y 

0  ->0 

1  2  -»Z 

x-^x 

C2(a) 

X  X 

y-4-y 

0  ->0 

1  z  -^z 

X-^T^+X 

C2(c) 

X  X 

d  -^4)  +  Tc 

y  -^y 

0  ->7I-  0 

z  z 

X  -^7t-  X 

1  C2(b) 

X  — >x 

(()  ->7r  +  (j) 

1 

i 

1 

0  ->7t-  0 

Z  Z 

X  -^-x 

In  particular,  we  find 

c  =  (-  l)  V^) 

c  =  (-  d'ijkM) 

The  relations  in  Equation  (10)  tell  us  that  vvc  cannot  use  the  expansion  (6).  Instead  we 
must  construct  a  basis  from  linear  combinations  of  J  ,±  K  ,M)  whose  symmetry  causes  them  to 

transform  the  same  way  under  the  operations  of  the  V  group.  If  we  define 

JKMs>  =  +  (-  K,M>]  for  K  0  and  s  =  0,1 

V  2 


(11) 


and 


|jOMs>  =  (-  1)'|J0M)  =  |jOM>forK  =  Oand  s  =  0 


(12) 


then  we  get 


EiJKMs>  =  |JKMs> 

Cj(J;KMs>=(-  1)''!JKMs) 

J-K+s 

C2(^JJKMs>  =  (-1)  IJKMs)  (13) 

C2(j,)|JKMs>  =  (-1)‘%KMs> 

Hence  the  factors  of  Equation  (8)  arc  divided  into  four  sub-blocks  and  the  eigenvectors  and 
eigenfunctions  of  each  sub-block  correspond  to  one  of  the  symmetry  operations  of  the  V  group  (see 
Table  5). 


Table  5 

Asymmetric  Rotor  Symmetry  Operations.  From  [25]  [p.  271], 


Secular 

Determinant 

Symmetry  Designation 

i 

1 

K  1 

s  i 

J  even 

J  odd 

E+ 

Even 

0 

A 

Bz 

E- 

Even 

1 

Bz 

A 

0+ 

Odd 

0 

Bx 

By 

0- 

Odd 

1 

By 

Bx 

Now  we  must  calculate  the  matrix  elements  of  that  appear  in  our  sub-blocks.  The 
commutation  rules  p  i'J  j]  “  ^ijk  ^  k  give  i 


!US 


(JKM|J  ,|J,K  ±  1,M>  =  i[  J{J  +  1)  -  K{K  +  D] 


OKMjJ  y|J,K  +  1,M>  =  t|(J(J  +  1)  -  K(K  +  l)] 


(14) 


<JKMIJ^JKM>  =  K 


which  imply 


<JKM|j^jK'M)  =  £<JKM!J  ,|JK"M><JK"Ki-J  ,jJK’M> 

K" 


(15) 


=  i[j(J  + 1)  -  +  |[J(J  +  1)  -  (K  ±  ])]  V(J  +  1)  -  (K  ±  1)(K  ±  2)]  ^8 


K,K'±2 


Similarly 


OKMJ^  JK'M>  =  -iDo  +  1)  -  K^]5 


KK- 


|[J(J  +  1)  -  (K  ±  D)  ( J(J  +  1)  -  (K  ±  1)(K  ±  2)]  ^5^  (16) 

lx  ,K  z2 


(JKM;j5;jK'M)  =  K“5 


and  finally,  for  an  asymmetric  rotor  in  the  prolate  limit 


KK' 


(17) 
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^KK'=  +  1)  -K^]+  AK2}5j 


+  [;j(B-C)[J(J  +  l)-K(K+:/) 


■•:j+l)-(K±l)(K±2)]  (5 


Note  that  for  a  prolate  symmetric  top  (B  =  C),  diagonal. 

Now  let  CJ  and  p  be  scalar  quantitk*.  lenwecansay 


Hr(oa*p,<,b*p.»c*p)  =  +  P)j'  +  p)];  +  (OC  +  p)J  ? 


=  a(AJ^,+  Bj^^  +  Cj;)  +p(j^,  +Jt -"Jc)  “oH 


R(  ABC) 


with  the  eigenenergies 


E/  „  _  ^  ~  u  \  Pj  (J 

(aA+p,0B+p,oC+p)  (abc) 


Now  let 


-(A+C) 
^  (A-C) 


which  imply 


'K,K’±2 
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oA  +  p  =  1 


We  begin  our  discussion  using  the  rigid  rotor-harmonic  oscillator  approximation.  This 
approximation  assumes  the  hamillonian  can  be  written  as 

(26) 

^The  theory  of  Section  2.2  is  taken  from  the  following  References.  [31  (pp.  94-95),  [25)  [p.  32);  (25) 
[p.82,2901. 
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and  therefore  the  vibrational  and  rotational  energies  as 


E  =  Ev+Er  (27) 

This  type  of  separation  also  leads  to  the  form  of  the  total  wave  function  as  a  superposition  of 
two  different  wave  functions: 


Vvr  =  VvVr  (28) 

To  find  selection  rules  and  line  strengths,  the  matrix  elements  of  the  electric  moments  connecting 
the  vibrational-rotational  states  V'R'  to  V"R"  within  the  same  electronic  states  are  considered. 
Microwave  studies  have  shown  that,  for  this  molecule,  the  permanent  electric  moment  lies  along 
the  axis  of  symmetry,  which  for  2  F  jup  representation,  can  be  designated  as  the  a  axis  (axis 
of  least  moment  of  inertia).  If  we  call  this  axis  the  z  axis,  then  both  sets  of  axes:  x,y,z  and 
X,Y,Z,  representing  axes  fixed  in  the  molecular  and  lab  frames  respectively,  are  needed  to 
formulate  the  problem. 

The  electric  moment  relative  to  the  body-fixed  z  axis  can  be  written  as 


Mz=(MJ  +2 


(29) 


Q^=0 


where  (M  ^  is  the  permanent  electric  moment  and  Q .  is  the  kth  normal  coordinate  of  the 
o 

vibrational  movement. 

The  component  of  the  electric  moment  along  the  lab-fixed  Z  axis  w'ould  be 


M  ^  =  M  X  COS  xZ  -r  M  ycos  yZ  +  M  j,  cos  zZ 


(30) 


V.  nere  the  direction  cosines  between  the  two  coordinate  systems  are  functions  of  the  rotational 
coordinates  alone.  Generalizing,  we  say 


99 


'^F  =  I‘^Fs(e,o,ri-'^8  ;  F  =  X,Y,Z  .  g  =  x,y,2 


(3i) 


O 


Fg(0,6,x)  the  direction  cosines  of  the  rotating  molecule  whose  molecule-fixed  frame  is 
described  by  the  Euler  angles  6i  <!)  >  X  • 

The  matrix  elements  of  the  electric  moment  are 


<VR'|Mp.V  ■«■•)=  |>i<V'?RMFVv.¥R.di:  (32) 

Substituting  Equation  (31)  and  the  generalized  form  of  Equation  (29)  into  Equation  (32)  enables 
the  separation  of  Equation  (32)  into  two  independent  parts: 


<V'R-:MpV'R">  =  X(Mg) 

g 


f33) 


and 


R  Fg 


;R"> 


(R* 

'^Pg  )  represents  the  matrix  element  in  the  F  direction  corresponding  to  a  unit  electric 

moment  along  the  g-direction  associated  with  the  change  in  rotational  quantum  numbers  from 
R '  to  R  " .  As  such  it  is  independent  of  the  vibrational  normal  coordinates  of  the  system. 
<v'*Mg;v")  are  the  matrix  elements  of  the  electric  moment  along  g  assodated  with  a  change 

in  vibrational  quantum  numbers  from  V  '  to  V  ” .  For  the  harmonic  oscillator  approximation  the 
eigenfunctions,  the  Henmite  polynomials,  arc  odd,  so  the  matrix  elements  are  only  non-zero  for  a 
change  in  vibrational  quantum  numbers  of  Av  ■=  ±  1 .  Such  transitions  are  called  fundamental 
transitions  because  they  are  the  most  intense.  To  obtain  overtone  transitions  for  which  Av  >  1 
or  combination  bands  in  w  hich  the  quantum  numbers  associated  v%ith  more  than  one  normal  mode 
of  vibration  changes,  it  is  necessary  to  include  higher  order  terms  in  the  normal  mode  expansion 
(29).  Coupling  between  normal  modes  caused  by  anharmonic  terms  is  zero  unless  v  =v"  (direction 
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cosines  are  even  functions).  Equation  (33)  represents  the  matrix  elements  for  pure  rotational 
spectra,  that  is,  transitions  between  rotational  states  in  the  same  vibrational  level.  These  types 
of  transitions  give  rise  to  microwave  spectra,  and  are  furictions  of  the  permanent  electric  moment 
of  the  molecule.Equation  (34)  represents  the  matrix  elements  for  simultaneous  changes  in  both 
vibrational  and  rotational  states. 

In  the  analysis  of  the  vibration-rotation  spectra,  absolute  intensities  are  not  necessary', 
only  the  relative  intensities  are  important.  Therefore  it  is  not  necessarv  to  know  the  exact  form 
of  <V'|M  ),  only  that  it  is  non-zero.  Hence  we  are  left  with  the  quantity 

2 

^|(R  jOpgjR  )j  ^  the  line  strength,  which  along  with  weight  factors  to  be  discussed  later, 
defines  the  relative  transition  intensities. 

2.3  Symmetric  Rotor  Selection  Rules^ 


Before  we  attack  the  transition  intensities,  we  first  address  the  selection  rules  by 
examining  the  matrix  elements  Pg  )  for  an  asymmetric  rotor.  The  selection  rules  for  a 

slightly  asymmetric  vibrating-rotating  molecule  can  be  best  understood  by  first  considenng  the 
selection  rules  for  a  symmetric  rotor.  This  is  because  for  a  s\  mmetric  rotor,  K  is  a  good  quantum 
number  (see  Appendix  B)  and  we  can  derive  exact  selection  rules  for  K  in  the  symmetne  limit. 

We  start  by  realizing  that  angular  momentum  will  transform  from  the  lab  frame  to  the 
fixed  molecule  frame  exactly  as  the  coordinates  transformed  under  rotation 


P 

Fg  (0,0.x)  F 


(35) 


Note  that  everything  is  done  in  the  representation  for  which  the  energy  of  the  symmetric  rotor 
is  diagonal  in  |  and  K.  In  this  representation  we  have 


^The  theory  of  Section  23  is  taken  from  Reference  (3)  {pp.  95-98); 


(36) 


OKM|PjrK'M')  =  K6„.5KK.5„„, 
and 

Okm|p1j'K-M’>  =  J(J  -f 

Consider 

[P^[P“,<I>]  ]  =  P^«D  -  2P^<DP‘  -f  OP'*  <38) 

Repeated  application  of  the  commutation  rules  pP  i'^j]  ~  ^ijk^k 


P’O  -  2P^OP^  +  OP"*  =  2h^(  P^O  -r  <DP^)  -  4h~P(P  ■  O)  <39) 

This  means  that 


(J’K’  M’Ip^O  -  2P^<D 


cDP%KM>=  2h^<J'K'MT^<D  +  OP^JKM) 


-  4h^(J’  K‘  M  ■  P(P  -  0)|JKM; 


(40) 


which  implies 


h1j’^I'+  1)'  -2J'(J-+  1)J(J  -i-  1)  -r  {-(J  OJKM) 


=  2h"[J'(J'+  1)  +1(1  +  1)]{J'K’M'  <1>  JKM>  -4h“a’K*M'P(P  -  0)JKM>  (41) 

For  J  the  last  term  on  the  R.H.S.  of  Equation  (41)  vanishes.  For  the  other  factor  on  the 
R.H.S.  of  Equation  (41)  we  write 

2[ J  (J’+  1)  +  J(J  +  1)J  =  (J'+  J  + 1)"  +  O’-  J)’  - 1  <«> 

The  term  on  the  LH5.  is  factored  to  give 
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•+ 1)^- 2j'0'-i-  i)j(j  - 1)  -i- + i)T =[(j'-  })\r* j - 1)] 


resulting  in 


[(J- +  J)  -i]0'K'M-O:JKM>  =  0  <«) 


or 


[(J’-r  J  -M)“  -  lj[o'- 1  f  -  1  _k J’K'  MloljKM>=  0  (45) 

For  a  non-vanishing  matrix  element  one  of  the  coetficient  expressions  must  be  zero.  Tiu;  first 
term  cannot  be  zero  because,  by  assumption,  JV  J  and  j’,]  ^  0 .  TIk:  sccorKl  term  is  zero  only  if 

r-j=±i. 

For  J  =  J',  note  that  P  -  <|)  is  a  scalar  otKjrator  and  is  therefore  invariant  under 

2 

operations  generated  by  P2aridP“.  Therefore  4vc  have 


[P;,,(p-<l>)3=[p^(p-o)]=o 

leading  to  the  following  matrix  representations 


(46) 


0  =  (J  K  M'£pI(P  -  <i>)]jKM)  =  (J  K  M-  ?;,(P  O)  -  (P  -  <D)P^JKM> 


=  (K*-  K)(J'K’M*lP  -  <!>JJKM)  (47) 

arKl 


0  =  {J'K'M*[pI(P  -  0)1kM>  =  (f  K  M-  P‘(P  O)  -  (P  -  0)P^JKM) 

=  Ir(J  -  D-  J(J-  l)l(J  K  -M  P  OJkM^  <4S) 
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Equation  (47)  and  Equation  (48)  lead  to  nonvanishing  matrix  elements  only  for  K  —  K  and 
J'=  J  .  Finally  we  construct  operator  ±  ,  which  together  withP^  and  the 

commutation  relations  used  above,  show  us  that 


(J'K’M’ 


d)  +  i(3>  ]KM)  =  55  5 

Fy  “  Fx  ^  M’,M±1 


(49) 


which  implies  the  selection  rule  AM  =  ±  1 .  Collecting  the  important  results  we  have  for 
the  selection  rules  for  the  symmetric  rotor  in  the  general  case  (V'  |M  g|V  ")  are 


AJ  =  0,±  1  ,  AK  =  0,±  1 


(50) 


2.4  Asymmetric  Rotor  Symmetry  Considerations  and  Selection 
Rules'^ 


The  spectrum  of  the  asymmetric  rotor  will  be  complicated  because  the  selection  rules  are 
more  complex  than  in  the  svTnmetric  case.  This  is  because  of  the  increased  number  of  energy 
lev'els  and  the  arbitrary  direction  of  the  changing  dipole  moments.  Despite  these  problems, 
general  statements  can  be  made  concerning  selection  rules  for  asymmetric  rotors. 

First,  since  total  angular  momentum  is  always  conserved,  even  though  the  molecule  is 
now  asymmetric,  the  transformation  matrices  for  the  asymmetric  case  are  diagonal  with  respect 
to  J  (as  they  are  in  the  symmetric  case).  Tlie  selection  rules  for  J  are  therefore  the  same;  that  is, 
AJ  =  0,±  1.  The  selection  rules  for  and  can  be  obtained  by  examining  the  symmetry 

of  the  molecule  in  question. 

The  rotational  behavior  of  a  molecule  of  C  2^  symmetry  depends  on  the  ellipsoid  of 
inertia,  which  is  symmetric  with  respect  to  a  rotation  of  180°  about  the  principal  axis  even 
though  the  molecule  itself,  being  asymmetric,  may  not  be  symmetric  with  respect  to  such  a 

^The  theory  of  Section  2.4  is  taken  from  the  following  References;  (2)  (p.  93];  [3]  Ipp.  106-107];  [8] 
[pp.  214-216]. 
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rotation.  The  rotational  wave  function  Vjk  K  is  either  symmetric  or  antisymmetric  with 
respect  to  such  a  rotation. 

For  an  asymmetric  rotor  in  the  limiting  prolate  case,  symmetry  of  the  wave  function 

with  respect  to  a  180°  rotation  about  the  a  axis  (axis  of  least  moment  of  inertia)  which  is  the 

±iK  X 

axisof  symmetry,  depends  on  the  angle  X  as  in  e  .  Therefore  ^JK  is  symmetric 

when  K_.j  is  even  and  antisymmetric  when  K_^  is  odd.  Similarly,  in  the  limiting  oblate  case, 
^  JK  ^  ’S  symmetric  with  respect  to  rotation  about  the  axis  of  greatest  moment  of  inertia 
(c  axis)  when  is  even,  and  antisymmetric  when  is  odd.  Since  successive  rotations  of 
180°  about  all  throe  axes  bring  the  molecule  back  to  its  original  position,  the  symmetry  for  a 
rotation  about  the  b  axis  is  the  product  of  the  a  and  c  axis  symmetries. 

Now  if  the  changing  dipole  moment  lies  along  the  a  axis,  the  matrix  element,  on  which 
the  relative  transition  intensities  depend,  is  of  the  form 


-1  1 


(51) 

whore  K'^  and  K".j  represent  the  quantum  numbers  of  the  initial  and  final 

rotational  states,  |i.  is  the  dipole  moment  and  cos(Fa  )  is  the  cosine  of  the  angle  between  a  and 
some  space-fixed  axis.  Since  cos(Fa)  changes  sign  for  a  180°  rotation  about  the  c  axis  for 
example;  then  must  also  change  sign  for  the  matrix  element  to  be  non¬ 
zero.  Otherwise  itself  would  appear  to  change  sign,  and  since  the  ellipsoid 

moment  of  inertia  is  symmetric  with  respect  to  this  rotation,  this  would  imply  the  matrix 
element  would  have  to  be  equal  to  zero  to  avoid  a  contradiction.  Therefore  transitions  of  this 
type  can  occur  only  when  K'^  and  are  of  different  parity.  A  similar  argument  for  rotation 
about  the  b  axis  shows  that  K'  and  must  be  of  the  same  parity.  This  type  of  procedure 

is  repeated  for  changing  moments  along  the  b  and  c  axes  to  establish  selection  rules  for  all  cases. 
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Symmetry  properties  for  asymmetric  wave  functions  are  summarized  in  Table  6  and  selection 
rules  in  Table  7. 


Table  6 


Symmetry  Properties  of  Asymmetric  Rotor  Wave  Functions.  From  [2]  [p.  94]. 


Behavior  with  180°  rotation  about  principle  axes 

a 

(least  moment) 

b 

(intermediate 

moment) 

c 

(greatest  moment) 

e  e 

+ 

+ 

+ 

e  0 

+ 

0  0 

- 

+ 

0  e 

- 

- 

+ 

Table  7 

Selection  Rules  For  Asymmetric  Tops.  From  |2)  [p.  94]. 


Axis  Parallel  to 

Changing 

Dipole  Moment 

—  _ :  __  r-:  - - *  *  »i - - 

Allowed  Transitions 

a  (least) 

ee  eo 

oo  oe 

b  (intermediate) 

ee  00 

eo  ^  oe 

c  (greatest) 

ee  oe 

00  eo 
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Briefly,  we  see  that  the  selection  rules  for  symmetric  rotors,  AK  =  0,±  1  have  been 
relaxed.  In  general  vve now  have  AK  =  0,±  2,±  4,...  or  AK  =  0,±  1,±  3,...  that  is,  AK_^ 
and  AK  j  are  restricted  to  even  or  odd  changes.  However,  not  all  numerical  combinations  of 
AK_.j  and  AK  ^  are  possible  since  K  _j  +  K  ^  =  J  for  even  levels  (even  levels  defined  as 
K_^  +  K^  +  J  =  even)  and  K _^  +  K .|  =  J  +  1  for  odd  levels  (odd  levels  defined  as 
K_^  +  K.j  +  J  =  odd). 

2.5  Relative  Intensities  of  Asymmetric  Rotor  Transitions  From  the 
Direction  Cosine  Line  Strengths^ 

2 

The  line  strengths  mentioned  earlier  were,  strictly  speaking, 

Fg‘ 

direction  cosine  matrices  of  basis  functions  for  a  symmctnc  rotor  described  by  the  Euler  angles.  To 
describe  the  line  strengths  for  an  asymmetric  rotor  we  just  represent  the  asymmetric  rotor  wave 
functions  as  a  sum  of  symmetric  rotor  basis  functions  which  transform  under  rotation  to  the 
degenerate  pairs  to  which  they  converge  to  in  the  limiting  symmetric  cases.  Hence,  direction 
cosine  matrices  for  the  asymmetric  rotor  can  be  calculated  from  linear  combinations  of  values  of 
d)  pg  which  give  the  same  symmetry. 

We  will  detine  transition  intensities  as  the  direction  cosine  line  strengths  as  modified 
by  the  fraction  of  molcv  '  available  in  the  ground  state  of  the  particular  transition.  In  the 
rigid  rotor-harmonic  oscillator  approximation  this  fraction  is  given  by 

'  ”  ^JK  K  (52) 

-1  1 


^The  theory  of  Section  2.5  is  taken  f  .n  the  foilowin..^  References;  12]  lpp.l(X)-101  ];  [3]  (pp.  105, 
110-m];[8][pp.  210-214]. 
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This  approximation  also  assumes  no  coupling  between  vibration  and  rotation,  which  is 
particularly  true  for  low  vibrational  quantum  numbers.  Note  we  have  also  assumed  all  the 
electrons  to  be  in  their  lowest  state.  The  fraction  of  molecules  in  a  particular  vibrational  state  of 
energy  hv(n  +  -|),  n  =  0,1,2,...  is 


exp 


fv  = 


-hvl 

V 

l) 

k 

T 

B 

Z  exp 
n=0 


-hv(n4] 


k  T 
B 


=  exp 


-n->pfs 


L 


(53) 


For  a  multi-band  harmonic  oscillator  it  is 


Ilexp 


m 


I"!  -  exp  ^ 


m 


A 


’^b’^  J 


(54) 


where  d  is  the  degeneracy  of  the  vibrational  mode  of  frequency  V  .  Since  there  is  a  2J  +  1 
degeneracy  from  the  M-quantum  numbers,  we  have  for  the  fraction  of  molecules  in  a  particular 
rotational  state  (the  asymmetry  of  the  molecule  splits  the  K-degeneracies  of  the  symmetric 
rotor) 


(2J+l)exp 

-E 

(abc) 

^B^  J 

X(2J+l)exp 

J 

’-E, 

(abc) 

J 

(55) 


where  E(abc)  is  the  energy  given  in  Equation  (25). 

In  our  case,  the  Vg  transition  is  non-degenerate  and  the  ground  state  has  the 

vibrational  quantum  number  n  =  0.  So  in  principle,  wc  can  write  the  intensity  of  a  spectral 
transition  for  a  single  band  in  the  rigid  rotor-harmonic  oscillator  approximation  as 
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^VR-,V"R"°' 


-'nv\  VR’(abc) 

ITT  fP  -- k-f—  2 

— — T — ^1^^'  R'|Mp|v’R">| 

V  n(abc) 

^gvR'^^P  k  T 
n  L  B  J 


where  V’R'  and  V"R"  stand  for  all  the  vibration.il-rotational  quantum  numbers  of  the  lower  and 
upper  states  respectively.  ^  V'R'  is  the  multiplicity  of  the  lower  state,  which,  since  we  assume 
all  the  electrons  to  be  in  the  ground  state,  includes  only  the  nuclear  spin  statistics.  E  ^  is 

the  asymmetric  rigid  rotor  energy  from  Equation  (25),  N  is  the  number  of  molecules  per  cubic 
centimeter  and  V  is  the  transition  frequency. 

To  compute  relative  transition  intensities,  we  note  that  the  denominator  and  the 
constants  in  Equation  (56)  are  the  same  for  any  given  molecule.  Also,  since  the  frequency  from  one 
end  of  a  vibrational-rotational  transition  band  to  the  other  does  not  vary  appreciably. 


1-exp 


^  kgTJ  is  essentially  essentially  constant  for  a  given  band.  And,  as  pointed  out  earlier, 

the  purely  vibrational  part  of  (V'  R'‘Mp  V"R")  is  not  important  for  calculating  relative 

intensities  as  long  as  it  is  not  zero.  Therefore  the  relative  intensities  in  a  vibrational-rotational 
band  can  be  computed  to  a  high  degree  of  accuracy  with  knowledge  of  E  ^ ,  the  Mne 

strengths  and  the  nuclear  spin  statistics. 

We  finish  the  discussion  of  relative  transition  intensities  by  addressing  the  issue  of  the 
nuclear  spin  statistics. 


2.6  Nuclear  Spin  Statistics  of  SF2O2® 


Anytime  a  molecule  has  two  identical  nuclei  (same  isotope  of  the  same  element)  which 
share  the  same  electronic  environment  v\ithin  the  molecule,  the  nuclear  spins  will  affect  the 


^The  theory  of  Section  2.6  is  taken  from  Reference  12)  |pp.  102-103]. 
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symmetry  of  the  overall  molecular  wave  function.  The  overall  wave  function  can  be 
approximated  as  the  product 


=  (57) 

where  e ,  ¥  v,  Vr  and  ¥  j  represent  the  parts  of  the  wave  function  dependent  respectively 
on  electronic,  vibrational,  rotational,  and  spin  coordinates.  Naturally  the  behavior  of  ¥-j- 
with  respect  to  any  symmetry  operation  depends  on  the  behavior  of  each  of  the  four  parts.  The 
electronic  wave  function^  for  most  polyatomic  molecules  are  in  the  ground  state  (and  therefore 
S)mrimetric)  at  room  temperature,  and  the  ground  vibrational  state  for  SF2O2  is  also  symmetric  in 
the  harmonic  oscillator  approximation,  so  the  symmetry  of  the  total  wave  function  for  this 
transition  in  the  SF2O2  spcctm  is  determined  only  by  the  product  of  the  rotational  and  nuclear 
spin  functions.  We  can  then  write,  for  symmetry  purposes, 

'^T“'^jK_.jK^ (58) 

If,  in  a  molecule,  two  equivalent  nuclei  occur,  the  molecule  will  have  a  twofold  axis  or 
symmetry,  and  consideration  of  the  nuclear  spin  statistics  with  respect  to  an  interchange  of 
identical  nuclei  coordinates  about  this  axis  must  be  considered.  In  SF2O2  vve  have  two  pairs  of 
identical  nuclei,  but  since  the  oxygen  nuclei  have  zero  spin,  they  will  not  contribute  to  the  spin 
statistics.  Only  the  two  fluorine  nuclei,  each  of  which  has  a  spin  of  1/2  ,  will  affect  the  nuclear 
spin  states. 

The  nuclear  spin  function  can  be  either  symmetric  or  antisymmetric,  and  is  formed  by 
taking  linear  combinations  of  the  spin  functions  of  each  individual  nucleus,  that  is 
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\l/j  =  an,(l)a^.(2) 


m  =  m 


(symmetiic) 


\j/,  =  a  Jl)cj_.(2)  +  a„.(l)a  J2) 


m 


m 


(symmetric) 


(59) 


V,  =am(l)cj  (2)-<y  (l)aj2) 

*  mm"'  (antisymmetric) 

where  G  j^(l)  and  G  ^^(2)  are,  respectively,  the  spin  states  of  the  first  and  second  identical 

nuclei  with  projection  m  of  spin  1  on  a  space-fixed  axis.  As  usual,  m  takes  on  the  21+1  values 

1,1  —  1,1  —  2,...,— 1.  Hence  there  are  21+1  combinations  possible  for  the  first  equation  of 


(59)  and 


1^(21-!- 


1)  -(21  rl) 


possible  combinations  each  for  the  last  two  equations  of  (59).  This 


gives  us  a  total  possible  number  of  symmetric  spin  functions  of 


N3yn,  =  (2I  +  l)(l  +  l) 

and  a  total  possible  number  of  antisymmetric  spin  functions  of 


(60) 


(61) 


A  molecular  rotation  of  180°  about  the  symmetry  axis  (a  axis)  in  SF2O2  amounts  to  a 
complete  exchange  of  molecular  coordinates  of  the  tw'O  fluorine  nuclei.  Since  only  the  fluorine 
nuclei  (which  are  fermions)  contribute  to  the  spin  statistics,  V  j  must  obey  Fermi-Dirac 
statistics  for  the  exchange.  Hence  we  must  always  pair  an  antisymmetric  nuclear  spin  function 
with  a  symmetric  rotational  function  to  keep  the  nature  of  ¥  j  antisymmetric.  Similarly,  a 
symmetric  nuclear  spin  function  will  always  be  paired  with  an  antisymmetric  rotational  wave 
function  for  rotation  with  respect  to  the  symmetry'  axis. 

Since  I  =  2  (or  the  fluorine  nucleus ,  Equations  (60)  and  (61)  tell  us  we  have  a  total  of 

three  symmetric  and  one  antisymmetric  nuclear  spin  functions.  So,  according  to  Equation  (56), 
rotational  wave  functions  for  SF2O2  marked  with  a  under  the  least  moment  column  in  Table  6 


34 


will  have  (all  other  factors  aside)  three  times  the  intensity  of  those  marked  with  a  These 
results  are  tabulated  in  Table  8. 


Table  8 

Nuclear  Spin  Statistics  for  SF2O2 


^-1  ‘S 

Symmetry  with  respect 

to  rotation  about 

axis  of  symmetry  (a  axis) 

Statistical  Weight 

e  5 

1 

i  e  0 

symmetric 

1 

0  0 

antisymmetric 

3  ! 

1  0  e 

antisymmetric 

3  i 
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Chapter  3 

Experimental  Procedure 


In  this  experiment,  a  tunable  diode  laser  was  used  to  generate  laser  light  in  the  11  ji. 
range  to  observ'e  the  SF2O2  spectral  profile.  The  laser  light  was  transmitted  through  a  150  cm 
absorption  cell  and  the  resulting  signal  was  chopped  at  500  Hz  before  entering  the  detector.  The 
spectral  trace  was  recorded  on  a  personal  computer  (PC). 


3.1  Experimental  Description 

The  layout  of  the  experiment  is  depicted  in  Figure  4. 
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Figure  4— Experimental  Layout 


The  source  is  a  small  (400  p.  x200  p.)  crystal  chip  of  the  Pb-salt  family  with  a  typical  line  width 

of  a  few  megahertz  (161  (p-  945),  (15)  Ip.2695l.  The  chip  is  mounted  in  a  refrigerated  dewar  (cold 
head)  and  is  part  of  an  overall  laser  system.  This  system  is  a  Laser  Photonics  model  LS-3  laser 
spectrometer.  In  addition  to  the  laser  source,  the  spectrometer  includes,  a  HgCdTe  type,  liquid 
nitrogen  cooled  detector  with  low-noise  amplifier,  a  cryogenic  temperature  stabilizer  (CHS  );  a 
laser  control  module  (  LCM  );  a  grating  monochromator;  and  various  optics  (see  Figure  5 ). 
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The  laser  is  cooled  by  a  CTl-Cryogenics  model  22C  cryodyne  refrigeration  ^slem.  It 
operates  on  a  high  purity,  closed  He  gas  expansion  cycle,  and  can  cool  the  iaser  down  to  a 
temperature  of  14°K  -  The  CTS  provides  temperature  stabilization  of  the  laser  cold  head  to  the 
order  of  less  than  a  millikelvin.  At  cryogenic  temperatures  it  also  provides  a  convenient  and 
accurate  temperature  indicator  for  the  range  of  operation  of  this  expenment.  Although  '.ne  cold 
head  is  a  closed  system,  interior  outgassing  necessitated  periodic  evacuation.  The  v'acuum  was 
maintained  by  a  liquid  nitrogen  cold  trapped  Varicn  diffusion  pump  operating  in  tande  ..  with  a 
standard  roughing  pump. 

The  LCM  provides  cur-ent  control  of  either  positive  or  negative  {xilarity  to  the  laser. 
The  current  can  be  adjusted  from  zero  to  a  limit  of  2  .Amps  at  a  continuously  adj  istable  rate  from 
10'^  to  10  Amps/sec. 

The  grating  monochromator  provides  two  functions.  It  gives  absolute  frequency 
calibration  to  within  05  cm  and,  by  varying  the  slit  width,  it  can  act  as  a  variable  frequency 
bandpass  filter  (-05  cm  w'de ).  The  latter  function  is  nccessaiy  because  the  laser  emits 
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multiple  modes,  and  high  resolution  spoctroscopv  is  performed  with  on  y  one  riKKie  at  a  time  1161 
Ip.  948). 

This  method  of  using  the  monochromator  and  the  CTS  for  coars?  tuning  and  the  LCM  for 
fine  tunirg  allov.":  ’  limited  tuning  range  of  about  i  to  2  cm'  * .  The  tuning  behavior  results  from  tiu? 
dependence  of  mode  frequenaes  on  both  the  energy  bandgap  and  tne  refractive  iruicx  of  tte  crystal 
1161  [p.9471. 

The  SI2O2  sample  is  stored  at  low  pressure  in  a  150  cm  stainless  sled  cylinder  with  AR 
coated  ZnSe  windows.  The  purity  of  the  sample  is  not  knoivn. 

Scans  were  recorded  on  a  computer  with  the  help  of  the  software  pad»agc  that  c*  >e  with 
the  Stanford  Research  Systems  model  SR510  Lock-in  amplifier. 

3.2  Conduct  of  the  Experiment 

From  earlier  experiments  we  knetv  the  P  and  R-branches  of  the  Vg  vibrational  band 

would  be  in  the  8S5  cm  region  111  1.  {HI ,  so  wt:  ordered  laser  crystals  which  were  designed  to  be 
tunable  in  this  region.  Once  the  laser  cry  stals  were  mounted  in  the  cold  head,  all  modes  of  the 
laser  output  in  this  particular  frequency  range  wise  identified  by  lunj:^  IIki  mortochromator  tt>  the 
approximate  wave  number  and  using  the  CIS  to  adpsi  tNt  temperature  of  the  laser. 

The  laser  output  was  obscr\  cd  on  the  osciilosGjpe  while  the  LCNl  ivas  used  to  ramp  the 
current  over  the  entire  current  range.  Only  iht>se  .»5Cillaling  laser  nwdes  whose  transmission 
profiles  were  free  of  mode-heps  or  otejcr  irregulanlies  were  selected  .  \Vc  tl^n  repeated  this 
process  at  different  tcmfxsralures  until  v.v  had  a  cuilcci:on  of  acccpsalrie  IIX^dcs  for  a  particular 
frequency  setting  of  tlw?  monochromatcr.  For  another  frequerKv ,  the  nKinochromator  was  rcsst  a:^ 
the  process  repeated. 

The  procedure  for  actually  recording  spectra  involved  three  steps.  .At  3  fixed 
monochromator  setting  and  a  fixea  cola  hvoa  temperature,  th;  laser  current  cummt  was  ^owiv 
swept.  First,  we  placed  the  cell  containing  the  SF2C)2  sample  in  the  ham  path.  The  ceil  tvas 


then  removed  from  the  beam  path  and  the  current  sweep  was  repeated  in  the  exact  same  manner. 
This  transmission  scan  provided  a  record  of  the  pure  laser  power  output  profile.  An  air-spaced  Ge 
fabiy-perot  etalon  (with  a  free  spectral  range  of  .05  cm  '^)  was  placed  in  the  beam  path  for  the 
final  scan.  The  etalon  would  hopefully  provide  a  means  to  establish  a  relative  frequency  scale  for 
each  family  of  scans  [16]  [pp.  948-949),  [13]. 

The  scans  were  recorded  in  a  personal  computer  in  a  volts  versus  time  Cartesian  format. 
This  computer  was  hooked  directly  to  the  lock-in  amplifier.  The  available  software  gave  us  an 
option  of  selecting  a  sample  rate  varying  from  20  to  2  Hz.  A  standard  sampling  rate  of  2  Hz  was 
selected  to  balance  the  need  for  a  high  signal  to  noise  ratio  against  the  need  for  a  reasonable 
scanning  time. 

Initially  the  SF2O2  sample  was  placed  in  a  cell  30  cm  long  at  300  mTorr.  This  cell  was 
small  and  easily  handled,  but  did  not  provide  for  enough  absorption  to  measure  low  J  rotational 
transitions.  The  SF2O2  sample  was  then  moved  to  a  cell  150  cm  long  and  kept  at  900  mTorr.  The 
pressure  was  gradually  reduced  to  150  mTorr  to  prevent  saturation  of  the  transition,  and  to  ensure 
we  were  well  within  the  doppler-broadenod  regime.  This  combination  of  cell  length  and  sample 
pressure  allowed  the  measurement  of  oven  the  lowest  J  transitions,  as  long  as  the  lines  were  not 
obscured  by  other  transitions. 

There  are  basically  two  methods  for  recording  spectra  with  a  tunable  diode  laser  [15]  [p. 
2697],  [13].  One  method  is  by  rapidly  ramping  the  current  over  the  entire  current  range  and 
recording  the  output  on  a  multichannel  integrator.  Since  a  multichannel  integrator  was  not 
available,  we  relied  on  the  more  conventional  lock-in  method.  In  this  method  the  diode  current  is 
slowly  swept  as  the  output  is  traced  on  the  recorder  (which  in  this  case  is  the  computer  screen) 
utilizing  a  step  and  integration  collection  method  [13].  After  some  trials,  we  decided  that  a 
current  scanning  rate  of  10’^  amps/sec  was  ideal  for  the  type  of  scans  we  wanted  and  the  amount  of 
detail  needed.  Faster  scanning  rates  tended  to  smear  away  details  while  slower  rates  allowed 
excessive  thermal  and  vibrational  noise  to  obscure  the  spectra.  The  characteristics  of  the  laser 
also  had  an  impact  on  the  data  collection.  A  crystal  with  a  frequency  scanning  rate  over  about  480 
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MHz/milliamp  would  force  us  to  slow  the  current  scanning  rate  below  the  level  at  which 
unacceptable  amounts  of  noise  would  creep  in. 


3.3  Data  Analysis 


After  a  set  of  scans  for  a  particular  monochromator  setting  was  recorded,  a  normalized 
spectrum  was  then  created  by  dividing  the  absorption  scan  by  the  transmission  scan.  It  was 
evident  that  the  scans  produced  a  distinct  set  of  lines  that  appeared  to  be  somewhat  regularly 
spaced,  but  before  those  lines  could  be  unambiguously  assigned  to  various  SF2O2  transitions,  it  was 
necessary  to  examine  the  available  information  on  the  molecule  more  closely. 

Previous  work  in  both  the  microwave  and  infrared  regions  has  given  reliable  information 
about  the  structure  and  inertial  parameters  of  SF2O2.  SF2O2  possesses  C2v  symmetry  with  a  pair 
of  oxygen  atoms  and  a  pair  of  fluorine  atoms  attached  to  the  central  sulfur  atom  tetrahedrally  (see 
Figure  6)  (14)  (p.  4).  From  previous  microwave  work,  the  three  inertial  constants  of  SF2O2  were 
found  to  be  A  =  .1713  cm"^,  B  =  .1693  cm'^,  and  C  =  .1686cm'^  (14).  The  shape  of  the  molecule 
suggests  (and  microwave  measurements  confirm)  the  assignment  of  the  least  axis  of  inertia  (a  axis) 
to  the  symmetry  axis  (z  axis)  of  the  molecule.  However  since  C  and  B  are  nearly  equal,  the 
assignment  of  the  b  and  c  axes  is  not  immediately  obvious.  Since  the  difference  between  the 
greatest  and  least  inertial  constants  is  small  compared  to  the  value  of  the  constants  themselves, 
SF2O2  is  classified  as  a  slightly  aspherical  asymmetric  top,  an  uncommon  type  of  molecule. 
Consequently,  the  asymmetric  splitting  will  be  small,  and  the  spectrum,  while  still 
representative  of  an  asymmetric  rotor,  will  be  simplified  114)  [p.  1). 

The  nine  normal  vibrations  for  SF2O2  are  given  in  Figure  7  along  with  their  symmetries. 
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Figure  7-The  Normal  Vibrational  Modes  Of  SF2O2.  From  (11]  [p.l792]. 
Figure  6  (previous  pagc)--Structurc  of  SF202-  From  (14). 


The  vibration  we  are  studying  has  been  described  as  the  SF2  asymmetric  stretch  ( V  g  band).  The 
schematic  of  the  transition  is  given  in  Figure  8. 


43 
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Figure  8--Schematic  of  the  Vg  Transition  for  the  Slightly  Aspherical  Rotor  SF2Q2 


The  spectrum  is  divided  into  P,  Q,  and  R  -branches  based  on  the  primary  selection  rules  for 
asymmetric  rotors  A  J  =  -1, 0,  +1  respectively  .  Because  SF2O2  is  very  nearly  a  spherical  top,  we 
should  expect  to  see  a  distinct  and  well-defined  PQR  structure  (11)  [p.  17911. 
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Classically,  a  vibrational  transition  is  connected  with  a  change  in  the  dipole  monaent  of 
the  molecule.  For  a  transition  to  be  IR  active,  it  must  have  non-zero  matrix  elements.  For  any 
molecule  of  the  form  X2YZ2  having  C2v  symmetry,  the  only  transition  that  is  IR  inactive  is  the 
tortion  oscillation  of  X2  against  Y2  (26]  [p.  240).  For  SF2O2  this  corresponds  to  the  mode  (see 

Figure  7). 

The  vibrational  transitions  of  molecules  are  classified  into  A,  B,  or  C  type  bands 
depending  on  whether  the  change  in  the  dipole  moment  lies  along  the  a,  b,  or  c  axes.  SF2O2  does 

have  a  permanent  dipole  moment  of  .228  ±  .004  Debye  which  lies  along  the  symmetry  axis  (a 
axis)  and  points  in  the  direction  of  the  two  fluorine  atoms  (14]  [p.  5).  Since  the  v .  asymmetrical 

stretch  of  SF2O2  lies  entirely  in  one  plane,  the  only  type  of  transitions  we  will  observe  will  be 

transitions  represented  by  changes  in  the  dipole  moment  along  the  a  axis  and/or  the  perpendicular 
axis  lying  in  the  plane  of  the  v.  vibration.  But  further  examination  of  the  molecular  structure  of 

SF2O2  suggests  that,  for  the  v  ^  asymmetrical  stretch,  the  strongest  dipole  moment  change  would 

lie  in  the  axis  perpendicular  to  the  a  axis  in  the  plane  of  the  vibration.  Since  earlier  microwave 
work  was  unable  to  resolve  the  b  axis  from  the  c  axis,  it  was  not  immediately  clear  if  we  would 
observe  B-type  or  C-type  bands  in  the  spectra.  In  any  case,  we  can  eliminate  A-type  bands  because 
the  structure  of  the  spectra  would  be  completely  different  from  what  we  observe.  Fortunately,  it  is 
relatively  easy  to  distinguish  between  B  and  C  type  bands  by  examining  their  Q-branch  structure. 
The  Q-branch  structure  of  a  B-type  band  docs  not  possess  a  strong  central  branch.  Rather,  the  Q- 
branch  has  a  central  minimum  and  its  lines  tend  to  overlap  the  P  and  R-branches.  In  contrast,  a  C- 
type  band  of  an  asymmetric  rotor  in  the  prolate  symmetric  limit  strongly  resembles  the 
perpendicular  band  of  the  symmetric  rotor.  This  band  is  characterized  by  a  strong  central  Q- 
branch  of  unresolved  lines  that  are  generally  stronger  than  tlie  P  and  R-branch  lines  [27]  [p.88]. 
This  stark  difference  in  Q-branch  structure  between  B  and  C  type  bands  should  help  to  determine  if 
the  observed  spectral  transitions  are  B  or  C  type. 

Finally,  we  can  use  Equation  (25)  to  estimate  the  relative  spacings  of  the  transitions,  at 
least  for  low  values  of  J.  If  we  neglect  centrifugal  distortion,  set  the  inertial  constants  for  the 


excited  state  approximately  equal  to  their  value  in  the  ground  state,  and  note  that  A  =  C,  we  see 
that  each  transition  is  separated  by  an  interval  E  =  (A  +  C)  =  10.2  GHz. 

So,  based  on  the  physical  properties  of  SF2O2,  we  are  looking  for  a  well-defined  PQR 
structure,  and  a  Q-branch  reflecting  the  characteristic  shape  of  a  (B  or  C)  type  transition.  The 
transitions  should  be  approximately  10.2  GHz  apart  from  one  another  and  each  transition  should 
have  a  fairly  regular  structure  due  to  the  slightly  aspherical  nature  of  the  molecule. 

3.4  Determination  of  Relative  and  Absolute  Frequency 


To  measure  the  spacing  of  the  transitions,  we  initially  used  an  air  spaced  Ge  etalon  with 

a  free  spectral  range  of  .05  cm'b  It  was  hoped  this  etalon  would  give  a  large  number  of  closely 

spaced  peaks  within  the  range  of  the  scan,  thereby  allowing  a  frequency  scale  to  be 

established[13].  Assuming  each  peak  was  separated  by  .05  em'b  a  linear  scanning  rate  for  each 

peak  interval  was  established  by  dividing  .05  cn.‘^  by  the  measured  time  interval  between  each 

etalon  peak.  The  spacing  of  each  transition  was  then  determined  by  counting  the  number  of  etalon 

peaks  between  the  transitions,  then  multiplying  the  determined  scanning  rate  by  the  remaining 

fractional  intervals.  We  were  able  to  find  that  the  interval  between  the  major  prominent  features 

of  each  scan  to  be  around  10  GHz,  leading  us  to  conclude  that  vve  were  in  fact  observing  transitions 
of  the  V  bandofSF202.  However,  we  also  measured  intervals  of  other  scans  with  the  etalon 

O 

that  ranged  from  8  to  12  GHz.  Apparently  the  free  spectral  range  of  the  etalon  was  too  large  to 
effectively  linearize  the  diode  scan.  Other  ctalons  with  smaller  free  spectral  ranges  were  tried, 
but  were  found  to  be  too  sensitive  to  vibrations  to  give  consistent  results. 

This  problem  was  corrected  at  the  same  time  a  way  was  found  to  determine  absolute 
frequency.  Since  the  monochromator  could  not  give  a  precise  location  in  frequency  space,  S1H4, 
which  has  a  structure  that  is  very  precisely  known,  vvas  used  as  a  reference  (28).  In  the 
experimental  procedure,  the  etalon  scan  v\as  replaced  with  a  simultaneous  scan  of  both  the  SF2O2 
sample  cell  and  the  SiH4  sample  cell.  The  SF2O2  transitions  \'-erc  now  superimposed  on  SiH4 
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lines  which  had  well  known  locations.  Relative  and  absolute  frequency  could  now  be  determined 
by  measuring  the  time  interval  between  the  nearest  two  SiH4  lines  which  straddled  each  SF2O2 
transition.  This  technique  did  not  completely  account  for  all  the  error  of  a  non-linear  diode  tuning 
rate,  but  result*?  were  much  more  self-consistent,  and  more  consistent  with  the  theoretical 
approximations,  (see  Table  14). 

3.5  Interpreting  the  Spectra 


An  examination  of  the  Q-branch  scans  showed  quite  clearly  that  we  were  dealing  with  C- 
type  transitions.  The  Q-branch  strongly  resembles  the  parallel  type  band  of  a  symmetric  prolate 
rotor.  This  knowledge,  along  with  the  nuclear  spin  statistics,  allowed  us  to  immediately  identify 
each  transition  wth  an  even  or  odd  rotational  quantum  number.  This  is  because  the  lines  of  each 
transition  of  the  P  and  R-branches  (the  Q-branch  lines  are  not  resolved)  at  some  point  show  a 
regular  alteration  in  intensity  in  complete  agreement  with  the  nuclear  spin  statistics  of  the  two 
fluorine  atoms.  For  a  C-type  band,  even  J  transitions  start  with  a  weak  line  (on  the  high  energy 
side),  odd  J  transitions  with  a  strong  line.  For  the  R-branch  these  parity  assignments  were 
confirmed  when  actual  J-values  were  assigned  to  each  transition  by  using  the  SiH4  lines  and 
counting  the  number  of  10.2  GHz  intervals  from  the  edge  of  the  Q-branch.  This  was  not  so  easy  on 
the  P-Branch  side  because  there  was  no  definite  edge  to  the  Q-branch  there,  and  it  was  clear  that 
an  unknown  number  of  P-branch  transitions  were  obscured  by  the  Q-branch  lines. 

To  be  absolutely  sure  of  the  quantum  number  assignments  on  the  R-branch  side,  and  to 
determine  how  many  P-branch  transitions  were  actually  hidden  by  Q-branch  lines  it  was 
necessary  to  actually  identify  in  the  scans  as  many  low  J  transitions  as  possible.  This  was  not  easy, 
because  the  low  J  transitions  on  both  the  P  and  R-branch  sides  were  weak  and  it  was  hard  to  tell 
them  apart  from  the  noise  in  the  scan.  One  way  to  strengthen  the  low  j  transitions  was  to  change 
the  Boltzman  distribution  of  the  sample.  We  did  this  by  cooling  the  cell  to  dry  ice  temperature. 
Under  the  approximations  given  in  Equation  (57)  and  Equation  (52),  the  rotational  wave  functions. 
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and  hence  the  rotational  energies,  are  independent  from  the  other  energies  of  the  molecule. 
Therefore  we  can  approximate  the  fraction  of  molecules  in  the  sample  at  a  particular  rotational 
energy  level  J  by  Equation  (55).  If  we  approximate  A  ~  C,  then  we  have  from  Equation  (55) 


(2J+l)exp  • 
X(2J+l)exp 

J 


-(A+C)J(J-fl) 


— - ^ — r  I  for  A  and  C  in  units  of  energy 

-(A+C)J(J+1) 


For  this  molecule  in  the  ground  state  »  ^(A  +  C) ,  even  at  dry  ice  temperatures.  Thus 
J  ,  and  the  denominator  of  Equation  (62)  becomes 


T  .  IN  -(A+C))()^f)  _ 

2kgT  (A+C) 


To  determine  the  preferred  J  value  at  any  particular  temperature  we  set 


which  leads  to 


J  max  ~  -J  2  for  A  and  C  in  units  of  energy  (61 

For  room  temperature  J  =  25-  For  dry  ice  temperatures  (-  80°C)  J  =20. 


By  cooling  the  cell  with  dry  ice  we  were  able  to  identify  all  R-branch  transitions,  and  all  P- 
branch  transitions  down  to  J=4.  (the  J  value  of  this  transition  was  later  confirmed  by  direct 
association  with  computer  generated  spectra.) 

Although  we  determined  the  SF2O2  spectra  to  be  primarily  C-type  transitions,  there  was 

still  some  structure  in  the  absorption  scans  that  could  not  be  accounted  for.  Since  there  were  no  B  or 

A  transitions  for  the  v  band,  this  additional  structure  could  only  be  coming  from  another 
8 

vibrational  transition  (hot  band),  or  from  an  isotopic  band.  The  only  isotope  that  exists  in  any 
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appreciable  quantity  for  any  SF2O2  component  is  Approximately  1/20  of  all  sulfur  is 
Therefore,  if  the  observed  structure  was  due  to  we  would  expect  to  see  a  C-t3q3e  spectrum  that 
was  1/20  as  intense  as  our  primary  spectra  with  the  following  modifications;  Under  the  harmonic 
oscillator  approximation,  the  vibrational  band  center  w'ould  be  shifted  by  an  amount  proportional 
to  the  ratios  of  the  reduced  masses;  that  is; 

_ _ 32  32 

ing+2m^+2mQ  ”  32+38+32  ~  102  (66) 


m , 


34 


ms+2fn  ^+201^ 


34+38+32 


34 

104 


(67) 


^%F2Q2  band  center:  (872.21)  » 869.12cm  ^  (5g) 

Furthermore  the  transition  spacings  will  be  shifted  an  amount  approximated  to  first  order  by 

,69) 

However,  this  tv-pe  of  structure  was  not  observed.  What  we  do  see  is  loo  strong  and  does  not  follow 
the  above  patterns.  Therefore  we  conclude  that  the  unidentified  structure  is  probably  not  due  to 
the  isotope 

To  check  the  possibility  that  the  structure  may  be  due  to  a  hot  band,  we  can  check  Table  9 
for  the  correct  vibrational  difference  combinations. 


49 


Table  9 

Observed  Frequencies  and  Band  Assignments  for  SF2O2.  From  [12]  [p.  572] 


I 

I 

We  look  for  difference  combinations  of  various  excited  bands  that  come  close  to  being  separated  by 
|[[  887cm'^,  and  have  an  appreciable  occupation  at  room  temperature.  The  only  likely  candidates  are 

listed  in  Table  10. 

I 

I 

I 

I 

I 

E 


Table  10 

Possible  Hot  Bands 


Band 

+  887  cm-1  = 

Band 

+  Left-over 
Rotational 
Energy 

Rotational  J- 
values^  to 
which  left¬ 
over 

Rotational 

Energy 

Corresponds 

V4=l 

(274cm-i) 

84cm"’  i  252 

'  (l077cm"’)  !  ! 

v.=  1 

D 

(360  cm“^) 

■  v,  =  l 

'  ( 1269cm-’) 

-  22cm"’ 

66 

V,=  l 
(386  cm-') 

12 

^  Based  on  a  spacing  of  approx.  .34  cm’^  per  transition 


Based  on  the  information  in  Table  10  we  can  eliminate  all  but  the  band.  The  other  modes 

involve  interactions  with  ]-vaIues  that  are  so  high  as  to  preclude  any  appreciable  occupation  at 
room  temperatureffor  example,  the  J=66  rotational  level  has  an  occupation  of  only  .028%). 
Therefore,  if  the  unidentified  structure  is  hot  band  activity,  it  is  probably  due  to  the  transition 

Next  we  use  the  fact  that  the  Boitzman  occupancy  of  the  band  will  change 
differently  with  temperature  than  the  occupancj  ef  the  Vg  band.  Using  Equation  (52),  Equation 
(54),  Equation  (62)  and  Equation  (63)  we  can  construct  the  following  table 
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Table  11 

Comparison  of  Fractional  Occupancies  (fyfj)  of  Typical  Low-J  Transitions  of  »he  Vg  Band  v/ith 
Hot  Band  =  1  Transitions  of  Similar  Er  argies 


Basic  Data 

Comparison  Of  Fractional 
OccupaiKies 

;  Vibrational 
j  State 

! 

Fractional  Fractional 

Occupancy  at  Orcupancy  at 

Room  Dry-  Ice  Temp. 

Temp.  (f  =  I93®K) 

(T  =  294®K) 

Fractional 
Occupancy  at 
l^m 

Temp. 

(T  =  294%C) 

Fractional 
Occupenq'  at 

Dry  Ice  Temp. 

(-T-  ^  I930g;) 

Ground  State 

fy  =  43.4% 

=  74.9% 

V,=  l 

=  6.6% 

t\.  =  4.2% 

.14% 

mmm 

iHlBnl 

Rotational 

State 

ground 

state 

.52% 

£  C  — 

grouna  j=/ 

state 

133% 

J  =  3 

Cj=.57% 

fp.87% 

^'^=1^=19  = 
.16% 

.13% 

J=7 

fj  =  1.19% 

fj  =  1.77% 

J=15 

fj=ZU% 

J  =  19 

mnn 

mm 

If  we  now  go  to  our  scans  we  sec  that  in  each  case  as  we  from  room  jcmperaiute  to  dry  ice 
temperature,  the  strength  of  the  v  barid  iran^tion  increases  while  that  of  the  hot  band 

O 

decreases,  just  as  pr«licted  in  Table  II.  If  we  lake  the  ratio  of  the  stron^t  v  tend  line  to  the 

S 


stron^t  hot  b-and  line  in  the  same  scan,  and  compare  this  rai,o  at  te  ih  Ujmra^turcs,  wv  should 
see  that  the  ratio  will  3h%'3ys  increase  as  we  go  from  rewm  icmpcralure  down  to  dry  ice 
temperature  Figure  9). 


the  C-type  transitions  of  the  Vg  band  are  most  probably  hot  band  transitions  originating  from 
the  Vg  vibrational  band. 
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Chapter  4 

ERROR  ANALYSIS 

The  largest  source  of  experimental  error  was  the  diode  itself.  Over  the  course  of  the 
scans,  taken  in  the  lock-in  mode,  small  deviations  in  frequency  output  linearity  were  observed 
which  did  not  reproduce  from  scan  to  scan  (15]  [p.2699].  These  instabilities  in  the  laser  affected  our 
ability  to  pinpoint  transition  frequencies,  that  is,  repeated  scans  of  the  same  transition  produced 
different  locations  in  frequency  space  for  the  same  lines.  The  precision  with  which  we  could 
measure  the  frequency  of  a  particular  transition  varied  from  transition  to  transition  depending  on 
the  behavior  of  the  laser  and  the  relation  of  the  transition  line  in  question  to  the  superimposed 
SiH4  sample  lines.  These  errors  wore  accounted  for  by  scanning  each  transition  several  times  and 
using  the  sample  standard  deviation  as  a  measure  of  the  uncertainty  or  precision  of  the  frequency 
of  each  line  entered  into  the  least-squares  fit  program.  The  program  then  assigned  a  relative 
weight  factor,  by  standard  least-squares  theory,  to  each  line  based  on  the  entered  uncertainty. 

The  instabilities  in  the  laser  output  frequency  arise  from  the  sensitivity  of  the  laser  to 
temperature,  current  and  mechanical  fluctuations.  In  each  case,  we  were  able  to  develop 
techniques  for  minimizing  the  errors  produced  from  each  source. 

Of  all  the  sources  of  instability,  the  diode  is  most  sensitive  to  temperature  fluctuation 
[15]  [p.  2695].  Temperature  fluctuations  o\er  the  course  of  a  scan  are  either  instantaneous  or  long 
term.  Long  term  temperature  fluctuations  (temperature  drift)  were  detected  by  comparing  the  CTS 
panel  readings  from  the  beginning  and  the  end  of  the  scans.  However,  unless  the  CTS  was 
monitored  constantly  during  the  scan,  only  cumulative  temperature  drift  could  be  detected.  For  the 
CTS  model  used,  which  employed  a  silicon  diode  temperature  sensor,  typical  temperature  drifts 
are  given  in  Table  12. 
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Table  12 

Temperature  Drift 


Temperature  Drift 

Si  Diode  Behavior 

T>25®K 

±.003-4^ 

mm 

2  7— 

Z../ 

T<25‘’K 

ol/ 

±.001^ 

min 

50  — 

Table  12  also  shows  that  the  temperature  is  susceptible  to  voltage  fluctuations,  particularly  at 
lower  temperatures.  Error  due  to  temperature  changes  was  minimized  by  rejecting  all  scans  for 
which  large  temperature  drifts  were  observed.  The  diode  temperature  was  also  allowed  a  period 
of  time  to  stabilize  to  the  new  current  setting  before  a  scan  was  initiated. 

Although  we  could  eliminate  a  large  number  of  scans  by  reading  temperature  drift  from 
the  CTS,  not  all  errors  from  the  temperature  fluctuation  could  be  eliminated  this  way. 
Instantaneous  temperature  variations,  temperature  hysteresis ,  and  temperature  drift  which 
affected  the  laser  but  could  not  be  measured  on  the  CTS  still  influenced  the  scans.  Temperature 
hysteresis  results  in  an  overall  (non-reproducablc)  frequency  shift  between  scans  [l.Sl  [p.  2699]. 

This  error  manifests  itself  when  the  scans  arc  normalized.  Since  absorption  and  transmission  scans 
were  not  taken  simultaneously,  dividing  them  together,  when  one  is  shifted  with  respect  to  the 
other,  may  have  introduced  a  shift  in  the  line  positions.  If  the  right  equipment  is  available,  the 
shift  between  absorption  and  transmission  scans  can  be  eliminated.  One  way  is  to  split  the  beam 
before  it  enters  the  sample  cell.  Part  of  the  beam  then  travels  through  the  sample,  while  the 
other  part  does  not.  The  beams  remain  spatially  separated  through  the  monochromator  and  enter 
separate  detectors.  Since  the  scans  arc  recorded  simultaneously,  there  w'ill  be  no  shift. 

But  even  if  absolute  shift  between  divided  scans  is  eliminated,  there  still  exists 
instantaneous  temperature  fluctuations  and  some  temperature  drift.  Both  drift  (which  has 
already  been  discussed)  and  fluctuations  are  functions  of  the  CTS.  The  ability  of  the  CTS  to 
stabilize  the  temperature  (and  hence  eliminate  fluctuations)  is  affected  in  part  by  the  sensor 
responsitivity.  The  increase  in  heat  capacity  of  the  sensor  with  temperature  will  slow  the  speed 
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of  the  cryogenic  heat  sink  and  increase  the  time  required  to  stabilize  the  temperature  as  the 
current  changes  during  the  scan.  We  compensated  for  this  as  much  as  we  could  by  optimizing  the 
gain  and  reset  controls  on  the  CTS  for  each  scan  temperature.  Another  way  to  eliminate  the  effect 
of  a  temperature  fluctuation  would  be  to  use  the  sweep  integration  method  of  recording  data.  Since 
scan  times  for  this  method  of  data  collection  are  extremely  fast  (10  seconds  or  less),  the  laser  mount 
temperature,  once  it  has  reached  equilibrium  with  respect  to  the  average  current,  generally 
remains  stable  over  periods  lasting  much  longer  than  the  typical  scan  time  [15]  [p.  2697]. 

Fluctuations  in  current,  or  non-linear  current  sweeps,  are  also  classified  as  instantaneous 
and  long-term.  Instantaneous  fluctuations  result  from  changes  in  ambient  temperature  and  the 
ability  of  the  LCM  to  control  transient  currents.  For  this  system,  typical  values  are  given  in 
Table  13. 


Table  13 

Current  Stability 


Current  Readout  Accuracy 

.ImAmp 

Drift 

mAmp  » 

=.02-r-^ 

day 

Temperature  Stability 

mAmp 

=.02-7^ 

Line  Voltage  Stability 

'  =  .01  mAmp  change  for  line  voltage 
'  change  up  to  20% 

From  the  table  w’e  see  that  the  current  is  much  more  stable  than  the  temperature.  But  at  a  typical 
diode  tuning  rate  of  480  MHz/mAmp,  even  current  fluctuations  can  lead  to  errors. 

The  major  source  of  mechanical  fluctuations  was  the  refrigerator  piston  in  the  cold  head. 
Although  various  methods  of  insulating  the  cold  head  were  tried  ,  this  vibration  could  not  be 
eliminated  completely.  This  is  one  of  the  reasons  the  etalon  w'as  abandoned  as  a  means  of 
measuring  relative  frequency  in  the  spectra.  The  jiggle  in  the  beam  as  it  traveled  through  the 
etalon  caused  the  peak  spacings  to  vary  considerably  within  the  same  scan,  making  them  almost 


useless  as  measuring  tools.  The  smaller  the  free  spectral  range  of  the  etalon,  the  worse  the  effect. 
Besides  going  to  the  SiH4  as  a  means  of  estimating  relative  frequency,  the  effect  of  the  vibration 
was  mininiized  by  using  irises  to  clip  the  perimeter  of  the  beam,  thereby  eliminating  the  most 
unstable  part  of  the  beam  profile. 

The  net  effect  of  these  error  sources  was  to  cause  the  current-temperature  combination  at 
each  frequency  to  be  different  from  scan  to  scan,  resulting  in  deviations  from  frequency  linearity 
down  to  the  MHz  level  (16]  [p.951].  Because  we  derived  relative  frequency  in  the  spectra  by 
estimating  a  linear  scan  rate  between  neighboring  SiHq  lines,  transition  lines  which  were  closer  to 
the  SiH4  lines  were  inherently  more  accurate  than  lines  which  were  farther  away.  Lines  that 
were  far  away  were  subject  to  a  cumulative  error  as  the  linear  estimate  deviated  farther  and 
farther  from  the  actual  non-linear  scan  rate.  This  effect  was  manifested  in  larger  sample 
variances  for  transitions  ’.hat  were  positioned  far  away  from  SiH4  lines. 

As  stated  earlier,  because  these  errors  manifest  themselves  as  variations  in  transition 
frequencies  which  are  not  reproducible  from  scan  to  scan,they  can  be  treated  as  random  errors. 
Therefore  the  sample  variance  is  a  measure  of  the  frequency  error.  Table  14  lists  the  data  from  the 
usable  scans  that  were  actually  entered  into  the  least-squares  fit  program.  It  was  assumed  that 
all  lines  within  the  same  transition  shared  the  same  variance.  Each  line  entered  into  the  program 
with  its  accompanying  variance  was  assigned  a  weight  factor  by  the  least-squares  routine.  Table 
14  does  not  represent  all  the  scans  taken.  Tliose  scans  which  were  too  noisy  or  otherwise  distorted 
or  smeared  beyond  value  were  eliminated  from  consideration  .  Scans  which  produced  values 
which  were  so  far  off  as  to  be  obviously  incorrect  were  discarded  or  retaken. 
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Table  14 

Spectral  Data  For  Least-Squares  Fitting  Routine 


Observed  J 

Transition 

Sample  Mean^ 

Sample  Standard 

Deviation 

Transition 

Spacing^  (Ghz) 

29  -4  30 

897.2524 

ti 

o 

o 

o 

>— a 

9.67 

28  -4  29 

896.9300 

.0005 

9.97 

27  -4  28 

896.5976 

.0010 

17^18 

893.2986 

.0002 

10.21 

16  -»17 

892.9582 

.0015 

12^13 

891.6268 

.0006 

10.13 

11-4  12 

891.2891 

.0004 

10.12 

10-4  11 

890.9517 

=.0001 

10.08 

9  ^10 

890.6156 

.0003 

10.09 

8  -49 

890.2792 

.0028 

10.20 
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Table  14  (cont) 

7-^8 

889.9391 

.0018 

10.20 

6  -^7 

889.5990 

.0016 

10.12 

5->6 

889.2617 

.0014 

10.18 

4  -^5 

888.9225 

.0005 

10.22 

3  -44 

888.5818 

.0006 

10.22 

2-43 

888.2412 

.0012 

10.21 

1  ^2 

887.9008 

.0001 

5  -44 

885.5030 

.0013 

10.20 

6  -45 

885.1631 

.0004 

10.23 

7-46 

884.8220 

.0005 

13  -4  12 

882.7302 

.0001 

16  ^15 

882.6873 

.0007 

^Reflects  value  of  first  line  on  high  energy  side  of  transition 

2  Rigid  Rotor  approximation  predicts  a  transition  interval  of  =  10.2  Ghz  for  low-J  transitions. 
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The  basis  of  accuracy  in  this  experiment  comes  from  the  SiH4  lines.  They  are  given  to  an 
accuracy  of  .0001  cm"^  [28].  Our  accuracy  may  have  also  been  affected  by  the  sample  rate.  Linear 
estimates  of  scan  rates  varied  typically  from  20-40  MHz/sec.  A  sample  rate  of  2  Hz  may  have 
introduced  errors  into  the  absolute  frequency  measurements.  However  the  net  effect  of  the  low 
sample  rate  over  our  entire  scan  should  cancel  itself  out,  and  this  effect  can  be  ignored. 

Since  we  w’ere  operating  with  a  sample  pressure  of  150  mTorr ,  we  were  well  into  the 
doppler-broadened  regime.  The  broadening  of  the  lines  is  given  by 


FWHM  = 


(70) 


where  FWHM  is  the  full  width  of  the  transition  line  at  half-maximum  in  MHz,  X  is  the 
frequency  of  the  laser  in  microns,  T  is  the  sample  temperature  in  °K,  and  M  is  the  mass  of  the 
sample  in  atomic  mass  units.  For  SF2O2,  v\  e  get  a  width  of  FWHM  =  57  MHz  at  room  temperature. 
This  amount  of  width  did  not  prevent  us  from  determining  line  centers  to  any  significant  degree. 

The  least  squares  fitting  routine  relies  on  the  ground  state  inertial  constants  as  starting 
points  in  the  fitting  process.  For  the  fitting  process,  the  ground  state  inertial  constants  from  the 
previous  miaowave  experiments  W'crc  used.  \Vc  can  use  the  rigid  rotor  approximation  and  our 
experimental  data  to  calculate  the  ground  slate  inertial  constants  and  verify  the  accuracy  of  the 
constants  used^.  We  start  with  the  general  rigid  rotor  hamiltonian  for  an  asymmetric  molecule 
in  the  prolate  limit 


which  means 


H.  =  +  Cj5 


(71) 


E(,bc)  =  'HR 


(72) 


^The  theory  of  this  portion  of  Chapter  4  is  taken  from  the  following  References:  [31  [pp.  43- 
44,188-189,204-205);  (29J  [p.  7381;  (91  [p.  15781- 
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where  {  )  stand  for  the  diagonal  elements  (average  values)  of  the  enclosed  operators  in  the 
basis  which  diagonalizes  • 

Now  let  H  be  perturbed  by  a  small  change  in  C,  denoted  by  5C  ,  and  let  the 
R 

perturbation  be  given  by  the  operator 


H’=5Cj5 


(73) 


Then  if  E  is  the  energy  eigenvalue  of  the  unperturbed  operator,  we  can  denote  the  first 

(30C; 

order  correction  due  to  the  perturbation  as 


E  =  (SCJc)  (74) 

^(a  b  c+Sc)  eigenvalue  of  the  hamiltonian 

H=  AJ^  +  BjJ-r(C  +  6C)J?  (75) 

and  define  an  e  such  that 


^(a/D.c+Sc)  ^(abc)^^  ‘  ^  (76) 

then 


But  since  e  =  0(5C^)  then 


°^(abc) 

8C 


lim 

6C-^ 


E'-j-e 

8C 


(77) 


dE 

tabc j  _,T 

3C  “ 


(78) 


In  this  vein,  we  can,  using  Equation  (25)  and  Equation  (5),dcfine  the  following  quantities 


62 


Table  15 

Approximate  Energ}'  Derivatives  of  a  Rigid  Asymmetric  Rotor  ifor  K~  -05).  '^rom  (3]  (p.  190J 


Ik  ,K,  !  “  :  P  j 

—  1  1  I 

Y 

330 

8.953 

1.838 

1.209 

440 

15.950 

2.368 

1.681 

550 

24.941 

2.936 

2.122 

551 

24.928 

2.990 

2.082 

541 

15.784 

8.668 

5.583 

542 

15.785 

8.492 

5.723 

532 

8.284 

15.048 

6.667 

533 

8.788 

12.297 

8.915 

523 

3.241 

20.880 

5.878 

524 

4.215 

11.508 

14.277 

514 

1.788 

16.959 

11.253 

515 

1.284 

4.713 

24.002 

505 

1.010 

5.452 

23538 

660 

35.932 

3.511 

2558 

661 

35.915 

3376 

2508 

651 

24.730 

10.321 

6.984 

652 

24.737 

10.280 

6.973 

642 

15.410 

16.817 

9.973 

643 

15.582 

15.990 

10.428 

633 

7.622 

25510 

8568 

634 

8.778 

18.881 

14341 

624 

3.328 

29388 

9.284 

625 

4.503 

15.434 

22.063 

615 

2.652 

20.147 

19.200 

6i6 

1.484 

5.830 

34.686 

606 

1.347 

6.219 

34.434 

770 

48.921 

4.087 

2.992 

If  we  no*w  use  the  correct  energy  difference  combinations,  we  can  isolate  the  ground  state 
constants.  Ground  state  energy  differences  are  given  by  the  difference  between  a  P-branch  and  an 


R-branch  transition  which  both  terminate  at  the  same  excited  stale  J 


1^  .  if  we  consider 


the  principle  sub-branch  of  a  C-type  transition  ,  then  each  transition  involves  ^  and 

AK  j  =  0 ,  and  the  ground  state  energy  difference  is  given  by 
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Col) 


where  AE  represents  the  energy  difference  between  two  ground  state  levels  appearing  in  the 
two  different  transitions. 

From  Equation  {25)  we  get 

AE=  (A  -s-C)(2J  -f- 1)  -f  J  (82) 

and  from  Equation  (SO)  ive  have 


+  [*^( K .K^)p  (86) 

Experimental  \-aIucs  of  AE  can  be  used  to  check  the  values  (to  ivithin  rigid  rotor  acturaoes)  of 
A,  B  and  C  determined  from  microwave  nrKssunsnents,  before  inputing  tlrem  into  the  least- 
squares  fitting  routine.  The  results  arc  tabulated  In  tint  rHJxi  chapter. 
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Chapter  5 


RESULTS  AND  RECOMMENDATIONS 


During  the  course  of  this  experiment,  we  scanned  over  an  interval  of  15cm"^  covering 
roughly  46  P  and  R-branch  transitions  and  thousands  of  lines.  The  200  lines  to  which  we  could 
•  eliably  assign  quantum  numbers  were  entered  into  the  least  squares  program  and  used  to  perform 
the  curve  fitting.  Over  4600  lines  have  been  mapped  using  the  S-Reduced  hamiltonian  to  model 
the  molecular  vibration-rotation.  Since  a  listing  of  all  4600  lines  would  he  prohibitively  long, 
graphical  representations  of  experimental  spectra  and  computer  generated  spectra  (in  the  P  and 

R-branches)  are  presented  for  visual  comparison. 

The  measured  excited  state  band  constants  and  V  g  band  center  are  presented  in 

Table  16. 


Table  16 

Measured  Excited  State  Band  Constants  and  Vg  Band  Center 


Band  Constant 

(em'h 

Calculated  value 

(cm'^) 

Estimated  Standard 

Deviation  (cm'^) 

''o 

887.2185 

.2513x10-3 

A 

1710 

.2376x10-5 

B 

.1691 

.2624x10-4 

C 

.1681 

.4281x10-4 

Dj 

.5809x10-7 

.6425x10-7 

D]K 

-.7969x10-7 

.9017x10-7 

dk 

.3439x10-7 

.3033x10-7 

Table  16  (cont.) 

dl 

.6288x10-7 

.1067x10-6 

d2 

-.2359x10-6 

.4440x10-7 

The  sextic  distortion  coefficients  were  also  determined,  but  were  found  to  be  so  small  as  to  not 
contribute  to  the  fit  to  a  significant  degree. 

Appropriate  energy  difference  combinations  of  the  most  reliable  experimental  low  J 
transitions  of  the  P  and  R-branches  were  used  to  check  the  values  of  the  ground  state  inertial 
constants  A,  B  and  C  determined  in  previous  microwave  experiments.  Equation  (85)  was  used  to 
calculate  A  and  C.  B  was  calculated  by  using  Equation  (86)  and  the  values  of  the  energy 
derivatives  from  Table  15.  The  results  are  presented  in  Table  17. 


Table  17 

Calculated  Ground  State  Inertial  Constants 


Inertial  Constant 

Calculated  Value 

(cm*^) 

Calculated 

Standard 

Deviation  (cm-^) 

Difference  From 

Value 

given  in  [14] 

(Mhz) 

A 

.17138 

.00018 

1.8 

B 

.16983 

.00080 

17.1 

C 

.16762 

.00221 

-28.5 

The  values  in  the  last  column  of  Table  17  are  a  reflection  of  the  fact  that  we  are  dealing 
exclusively  with  C-type  transitions  for  this  vibrational  band.  Because  for  the  principle  sub¬ 
branch  of  a  C*type  transition,  ^  ‘"^nd  ^  =  0  ,  vve  should  expect  the  transition  energy 

differences  to  be  primarily  deper  *  on  A  for  a  molecule  approaching  the  prolate  symmetric 
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limit  [see  Equation  (86)).  For  more  reliable  information  on  B  and  C  we  would  have  to  search  for 
those  C-type  transitions  for  which  AK  ^  >  0 .  These  are  the  weaker  transitions  in  a  C-type 

spectra  and  would  be  correspondingly  more  difficult  to  isolate.  That  is  why  the  error  differences 
for  B  and  C  in  Table  17  are  so  much  greater  than  the  error  difference  for  A,  and  the  estimated 
standard  deviations  in  Table  16  are  similarly  skewed.  Nevertheless,  the  differences  are  still 
within  our  experimental  resolution,  and  wo  can  use  the  values  of  the  ground  state  inertial 
constants  derived  in  Reference  [14]  in  the  fitting  routine  with  full  confidence. 

RECOMMENDATIONS 

The  TDL  is  still  a  practical  tool  for  high-resolution  spectroscopy  in  the  middi^  IR 
region.  However,  to  reduce  error,  the  sweep-integration  method  of  data  collection  is  preferable 
to  the  lock-in  method.  Since  the  data  points  are  integrated  simultaneously,  even  fast  scans 
would  be  long  enough  to  reduce  the  noise  and  prevent  line  smearing,  thus  giving  spectra  that  are 
reproducible.  Experimental  setups  using  the  multichannel  integrator  can  be  found  in  Reference 
[15]. 

TTie  next  logical  step  in  studying  SF2O2  should  involve  hetrodyning  a  CO2  laser  with 
the  diode  laser.  The  P(20)  line  of  ^^€02  (at  896.90947969  cm’bis  very  close  to  the  R-branch 
transition  J  =  28  — >  29  and  can  be  used  to  pinpoint  a  line  very  exactly.  This  type  of  experiment 
would  provide  a  very  reliable  method  to  independently  confirm  the  results  of  our  experiment. 
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The  following  8  pages  contain  graphical  output  of  experimental  spectra  (top  of  page) 
and  corresponding  computer  generated  spectra  (bottom  of  page).  The  transitions  depicted  here 
are  representative  of  the  P-branch  and  R-branch  transitions  that  were  scanned  during  the  course 
of  the  experiment.  For  the  sake  of  clarity  some  transition  quantum  numbers  and  frequencies  have 
been  ommited  in  the  experimental  spectra. 
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wavenumber(cm-1 ) 


v/avenumber(cni-1) 


absorption 


’-6 
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Appendix  A 


Development  of  the  Hamiltonian 

To  find  the  vibrational-rotational  energies  of  an  asymmetric  molecular  rotor  we  use  a 
model  which  treats  the  nuclei  as  point  masses  moving  through  a  potential  field  created  by  the 
average  motion  of  the  molecular  electrons.  The  molecules  are  thought  of  as  being  held  together 
by  semi-rigid  springs  whose  energies  arc  eigenvalues  of  the  Schroedinger  equation. 

HV  =  E¥  (A-l) 

A-1  The  Classical  Hamiltonian^ 

The  form  of  the  hamiltonian  is  found  by  first  considering  the  classical  kinetic  energy  for 
the  system  of  point  masses.  The  most  convenient  set  of  coordinates  to  use  consrats  of  a  ^t  of  axes 
fixed  in  space  ( X,Y^ )  and  a  set  of  axes  (  x,y^ )  rotating  with  the  molecule  whose  origin  is  at 
the  molecular  center  of  mass  and  whose  orientation  in  space  is  described  by  the  three  Eularian 
angles  ( 0,0,% )  -  The  Eularian  anglK  are  discussed  more  fully  in  Appendix  B. 

Let  the  position  of  the  nuclei  with  respect  to  the  molecular  center  of  mass  be 

I 

with  omponenfs  ( ^  a’  ^  a  ^  The  center  of  mass  is  located  at  a  distance  ir|  from  the  fixed 

axes.  The  equilibrium  position  of  the  particle,  which  is  fixed  in  the  moving  system,  is 
given  by  a  ^  (see  Figure  A-l). 


%he  theory  of  Salion  A-l  is  taken  from  the  following  References;  PI  (pp.  2-61;  HI  Ipp- 14-17, 
273-2781;  Wfpp.aW)-^!. 
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T  =  a(  V tot ) ^  =  yX +  (^  X  r  J  ■  ( H)  X  r  „)  +  V „  +  2  R  •  (cb  X  ?  J 

a  a 


-Z  ik- v„+2cb-  (r„x  vj]  (A-5) 

Since  vve  have  defined  tl.e  (  x,y,z  )  axes  to  be  a  center  of  mass  system,  it  is  true  by  definition  that 

Xm„r„=0  (A-6) 

a 

Taking  the  total  differential  of  is  a  function  of  rotation  and  vibration  )  with  respect  to 

time  gives 


0  =  ?^  =  2m  J(cb  X  r  J  +  V  J  =  5)  X  =  X^ccV^  =  0 

a  a  a  a  a 

(A-7) 

Since  the  atoms  in  the  molecule  are  all  vibrating  about  their  equilibrium  positions,  the 
rotating  axes  cannot  be  attached  to  the  rotating  coordinate  system  in  a  definite  way.  However, 
to  ensure  the  rotating  axes  actually  rotate  with  the  molecule  we  can  impose  the  condition 


Xmoc(a„  xrj  =  0  (A-8) 

a 

This  is  almost  the  .  me  as  stating  that  there  must  be  no  angular  momentum  with  respect  to  the 
body-fixed,  rotating  system  of  axes.  We  take  the  total  differential  of  Equation  (A-8)  with 
respeci  io  time: 


0  =  Xma[(a^x  r„  +  Va)  +(a„x  v„)]  =  2m„((b  xa„)  X  2mt,a„  x  (o)  x r„) 

a  a  oc 


a  a 


(A-9) 
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Using  Equation  (A-3),  Equation  (A-6),  Equation  (A-7),  and  Equation  (A-9)  to  modify  Equation 
(A-5)  leaves 


2T  =  i  M  +  2^m„(a)  xr„)- ((bxr„) +5^mcc(Va)  +2tb(paXVc)  (A-10) 

a  a 

where  M  =  ^ 

Disregarding  the  translational  term  and  noting  that  by  expanding  the  determinant  form  of  the 
cross  product  we  get 


Xm„(cbxr„) 

a 


(“Xra)  =  Ixx"x+Iyy®y+l2z“2-2Ixy“x"y-2I^C0vC0 


xy 


yz~y~z 


21 

and  the  kinetic  energy  becomes 


(A-11) 


2T  ~  I  JQ,0)  +  I  yyCO  y  +  I  b)2  2Ijg(0y(02  21  jjy  CO  0)  y  2Ij|2C0j(C02 


+  En^aVa  +  2C0x£m„(p„X  V„)  +  2c0yXma(p„  X  v„)  +  Zco^^niaCpa  x  Va) 

''  y  a 


a 


a 


o 


(A-12) 


where  Ixx/  lyy/  Izz/  ^xy/  lyz/  ^xz  ^re  the  components  of  the  inertia  tensor 


1  = 


XX‘ 


xy^  xz 


I  yxl  yyl  y2 
zx  ^  ^  zz  y 


(A-12a) 
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To  deal  with  the  coupled  terms  in  Equation  (A-12),  it  is  necessary  to  treat  the  vibrations 
of  the  atoms  as  smalj  displacements  about  their  equilibrium  points.  The  kinetic  energy  is  then 
given  by 


N 

S. 

a=l 


r  2 

2 

21 

(?“)  +1 

L\  dt  / 

<dtJ 

^  dt  y  J 

(A-13) 


where  Ax  ^ ,  etc.  are  the  components  of  P .  Let  us  now  replace  the  coordinates  Ax , , Az  , 

I  N 

by  a  set  of  mass  -  weighted  coordinates  q^  , ... ,  q^^^  defined  by  q  j  =  Ax  ^  ^ 

2  ~  a/^I  1 '  3  ~  a/^  1  ^  1  /  4  ~  a/^  2  ^^2 '  kinetic  energy  now 

becomes 


3N 

2T=Iq? 

i=l 


(A-14) 


The  potential  energy  V,  a  function  only  of  the  small  vibrational  displacements,  is  now  a  function 
of  the  q's.  We  can,  therefore,  expand  V  as  a  power  series  in  the  qj’s. 


2V=  2Vg  + 


3N  SN/  2  ^ 

q,q 

j_]V  lyq.  o  i,j=lV  *  j/q.=q  =0 


(A-15) 


Adjusting  the  zero  point  of  the  energy  to  be  the  equilibrium  point  and  noting  that,  at 


fdW 

equilibrium,!  9q . 


=  0 


,  v/ehave 


q.=0 


3N 


S^ijqjqj  ,  where  - 


i,j=l 


q.q.=  0 
>  ) 


In  deriving  Equation  {A-16),  the  higher  order  terms  have  been  neglected. 


(A-16) 
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Newton's  equations  of  motion  are  written 


d  ar  ,  av  » 

dt  dq.  +  dq.-^ 


,j=l,2,...3N 


(A-17) 


since  T  is  a  function  of  the  velocities  only,  and  V  is  a  function  of  the  coordinates  only. 
Substituting  T  and  V  from  Equation  (A-14)  and  Equation  (A-16)  into  Equation  (A-17)  gives 


3N 


qi+Sfiiqi=<>  ,H.2 . 3n 


)  M 

A  solution  to  this  set  of  3N  simultaneous  second-order  differential  equations  is 


(A-18) 


q.  =  A..cos|^X 


i  I 

+  £7 


(A-19) 


By  substituting  Equation  (A-19)  into  Equation  (A-18),  the  following  set  of  equations  result 


3N 

X(  f ij  -  5  j  A)  A  .  =  0  ^  j=l,2, ...  ,3N  ( a-20) 

where  ^  ij  is  the  Kronecker  delta.  Non-trivial  solutions  for  Aj  correspond  only  to  those  values 
A,  for  which  the  following  secular  equation  is  satisfied 


(^1,1 

^,3N 

T 

^2,1'(^2,2~  ^)'^2,3 

^2,3N 

f 

! 

=  0 

(A-21) 

^3N,1'^3N,2'^3N,3'" 

'(^3N,3N 

)! 

When  we  choose  a  certain  value  A.  =  X  so  as  to  make  the  above  determinant  vanish, 
the  coefficients  Aj  =Aii<  are  determined  for  that  particular  value  .  Since  the  above  system 
of  equations  docs  not  determine  Aji^^  uniquely,  but  gives  only  their  ratios,  arbitrary  solutions 
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(which  are  unique)  can  be  designated  by  the  quantities  Ij];  '"hich  are  defined  in  terms  of  an 
arbitrary  solution  A’jk  as  follows 


A* 


ik 


2 


(A-22) 


The  advantage  to  this  is  that  the  amplitudes  are  normalized  such  that 


Slik  =  1  (A-23) 

i 

Note  the  secular  equation  (A-21)  consists  of  3N  rows  and  3N  columns  corresponding  to  3N 

unknown  Aj's  and  3N  unkno^vn  eigenvalues  A.  •  There  are  however  six  roots  which  are  zero, 

k 

corresponding  to  the  three  normal  coordinates  of  fanslation,  and  three  normal  coordinates  of 
rotation,  leaving  only  3N-6  non-zero  roots  (corresponding  to  the  3N-6  normal  vibrational  modes). 
This  fact,  along  with  the  normalization  relation  (A-22),  allows  us  to  define 


3N-6 


m „  Ax ^  X^cck*^k 


k=l 


3N-6 

k=l 


(A-24) 


3N-6 


k=J 


k 


Jl 

V 

a  a 


=  Xa 

a 


+  A  yl  +  ^  'zl, 


2 

k 
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Now  the  coupling  terms  become 


a 


Xm  jAy„Az„ 

a 


3N-6 


I3,Q 


k 


Xm„(p„xvJ  = 
a  y 


3N-6 


k 


(A-25) 


Xm„(p„x  vj 


3N-6 

= 


k  ^k 


where  X^a(^ctl^ak  ^al‘ak)^l 

a,l 

^k~  2]r^a(^al^ak“^al^ak)^l  (A-26) 

a,l 

K|=Tm„fi  ,r  ,-r  ,1  .')Q, 

arwH  k  ^  “V  a!  ak  al  ak/  1 

^  a,l 


The  kinetic  energy  is  now 


^  ^  xx®"x  ^  y  ^  ^  zz^  z  ^^xy®x®y  ^^yz^y^z  ^^>z^x®z 

+  2(0  k  ^  k  vX^  k  Q  k  k  Q  k  Sq  k 

k  k  k  k 

We  will  find  it  convenient  to  have  the  kinetic  energy  expressed  in  terms  of  angular 

momenta,  rather  than  angular  velocities.  The  total  angular  momentum  is  a  vector  whose 

properties  arise  from  its  definition 
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?a=  X  («X  f  J]  +  Xm„(f„  X  T/J  (A-28) 

OL  a  a 

9T 

By  expanding  Equation  (A-28)  or  by  noting  P  ^  ,etc. ,  we  have 

^  X  ~  ^  XX®  X  “  ^  ""  ^  >z®z  X^  k  ^  k 

Py=-I^“x  +  I„»y-I)*“z+X5i^Q^  (A-29) 

Pz  =  -'>z®x-I)z“y+Izz“z  +  Z!'kQK 
The  momentum  conjugate  to  Qk  is 

Pk^"^  ~^k'^^k“x  ^k“z  (A-30) 

k 

Solving  Equation  (A-30)  for  ^  and  substituting  the  result  into  Equation  (A-29)  yields 

Px  “  I  XX®  X  ~  ^  jy®  y  “  ^  x2®2  X^  k(  ^k  “  ^k^’^  ~  ^  k^y  ~  ^  k^z)  (A-31) 
It  is  now  convenient  to  introduce  the  following  definitions 

^  xx~  ^  XX  “  X^  k  .  I  xy'^  ^  xy  X^  k^  k  •  P  ~  X^  kP  k 

^  yy“  ^ yy  “  X^ k  ;  ^  >z=  I >z ■*■  X^ k^  k ;  P y  “  kPk  (A'32) 

^  zz~  ^zz  “  X^  k  •  ^  ^  xz  X^  k^  k  •  P  z~  X^  kPk 

We  can  see  from  their  definition^  that  Py  p^are  functions  of  the  vibrations  alone.  As  a 

result  they  are  often  called  components  of  internal  or  vibrational  angular  momentum. 

Substituting  Equation  (A-32)  into  Equation  (A-31)  gives 
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(^x  P  x)  ^  xx^x  ^  xy^y  ^  xz®z 
(Py-Py)  =-I  I  ^0y-I  (A-33) 

(Pz  “Pz)  xz^x”^  jz^y”^!  zz 

The  inverse  transformation  of  Equation  (A-33)  is 


®x  Pxx(Px  Px)'^Pxy(Py  Py)‘^Pxz(Pz  P  z) 

“y  “  Pyx(P  X  ~  Px)  ^yy(^y  ~  P  y)  PyzC^z  “  P  z)  (A-34) 


^z  Pzx(Px  Px)'^Pw(Py  P  y)  p2z(Pz  “Pz) 


where 


.1  |■2  ji  I.  I.  . 

yy*  zz"‘  yz  ‘  zz*  yz*  xz 

Pxx=— -T-  ■Pxy=  - Z - 


Pxz 


i  Jiyl  yZ+1  yyl  5g, 


pyz 


^  xx^  xy^  >z 


Pw  = 


r  r  -!' 

xx‘  2Z  ‘  XZ 


I’  r.„-r 


.2 


.  Pzz  “ 


»:  >y  xy  (A-35) 


‘T*  —  T*  T* 

r  XX  y  ^  xy »  ”  ■  X^i 
I 

_ T»  r»  _  T*  • 

^  xy»^  yy»  ^  yz; 

-  I  ’  -  r  r 

^  XZ  »  *  \z  >  zz  • 


Substituting  Equation  {A-32),  Equation  (A-33),  and  Equation  (A-34)  into  Equation  (A-27)  shows 


the  kinetic  energy  as 
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2T  ^lxx(Px  Px)  ■^P;A(Py  P  y)  +Pzz(Pz-Pz)  +  xy(Px  “  P  x)(Py  “  P  y) 


j,  (A-36) 


+  2Hj2(^y  Py)(Pz  V 7x.iy z  Pz)(Px“  Px)  “XP^ 

Note  that  the  Pgg.  are  functions  of  the  normal  (vibrational)  coordinates  only.  The  classical 
hamiltonian  is  complete  when  V  is  added  to  Equation  (A-36). 


A-2  The  Quantum  Mechanical  Hamiltonians 


Because  the  total  angular  momenta  Px,  Py,  P2,  in  Equation.  (A-36)  were  not  derived 
with  respect  to  any  conjugal  coordinate  system,  we  cannot  transform  directly  from  classical  to 
quantum  mechanical  hamiltonian.  Thus  we  must  derive  the  quantum  mechanically  correct  form 
of  the  hamiltonian  in  an  arbitrary  coordinate  system. 

In  most  normalized  coordinate  systems,  an  element  of  differential  length  in  n  dimensions 
is  characterized  by  the  expression 

ds2=h,,dqf-rh22dq|-}--  ••+hn„dqj  (A-37) 

For  an  arbitrary  coordinate  system  of  n  dimensions,  a  differential  clement  of  length  is  given  by 


ds^  = 


2  du^ 

61 


2g,2du,du2  +  . 


^32^*“  2  ^8  23'^ '*2'*  “3 


onn^*'^n 


(A-38) 


where  the  coefficients  gjj  may  be  functions  of  the  coordinates.  The  transformation  of  the 


Schroedinger  equation  for  a  sing'e  particle  from  the  n-dimensional  coordinate  system 


^The  theory  of  Section  .4-2  is  taken  from  the  following  References;  (5j;  (1)  (pp  279-'’831‘  161 
[pp.262-2631. 
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represented  by  Equr ‘ion  (A-37) ,  to  the  n-dimensional  coordinate  system  given  in  Equation  (A-38) 


where 


yg2  3 

i-j  o  0u  o  o  3u , 

i.j=l  ‘V  1 


^^(E-V,)v,=  0 

n 


{A-39) 


JJ-  ■■  Jvq¥qg  ^  ciUjdu^  ■  •  -d  u  „  =  I  and 


The  corresponding  normalization  for  V  u  is 


|§1  l'Si2 

i§21»822 


!Snl’On2’  ’»gnn 


’Sin 

■S2n  (A-40) 


JJ  -- JVuV‘udu,du2-  -du„=l 


(A-41) 


Equation  (A-40)  and  Equation  (A-41)  imply 


VQ=g'’¥u 


Substituting  from  Equation  (A-42)  into  Equation  (A-39)  and  multiplying  through  by  _ 


(A-42) 

-1 

1,2  T 

h  g 

2m 


n  I  (  zi  ..  1  \ 

S^s'^Y^s'^Vu  +V„Vu-E¥u  =  0  (A-43) 

i,j=i  iV  j  y 

hA 

Replacing  the  differential  operators  ,  by  the  corresponding  momenta  and  absorbing  the  m 

into  the  determinant  gives  the  hamillonian  in  the  correct  form  fer  arbitrary  coordinates 
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{A-44) 


I^=lg'‘SPiS^g‘’Pjg'*  + V 

i.) 

Notice  that  the  potential  V  is  unchanged  from  its  classical  form  because  it  is  a  function  of  the 
coordinates  and  not  of  the  momenta. 

For  the  hamiltonian  in  Equation  (A-44)  the  kineHc  energy  is  given  by 

2T=Xgijqiqj  (a-45) 

'.j 

or  in  terms  of  the  momenta 

2T=i;g‘^PjPj  (A.46) 

i.j 

The  coefficients  g*J  are  elements  of  the  matrix  which  is  the  inverse  of  the  matrix  formed  by  the 
gjj,  and  the  momenta  pi  are,  of  course,  conjugate  to  the  generalized  coordinates  q,. 

Transforming  now  to  a  set  of  momenta  Pm  which  are  not  conjugate  to  any  set  of 
coordinates,  but  are  instead  denned  by  the  expression 


(A'47) 


now 


2T=  XG  “P’mP’ 


m,n 


{A-4S) 


where  we  define 


,mn 


=ys 

^  kmo  ^In 

k,l 


(A-49) 


We  must  now  determine  under  what  conditions  a  quantum  mechanical  hamiltonian  may  be 
written  in  a  form  analogous  to  equation  Equation  (A-44) ,  but  in  terms  of  the  momenta  P'nv 
namely 
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(A-50) 


1  zl  ± 

H=^G"2:F^G"'"G"r„G%V 

m,n 

where  G  =  !g  I . 

To  find  these  conditions,  substitute  Equation  (A-49)  into  Equation  (A-50)  and  the  expression 

P'in=Es"'‘P5  (A-51) 

i 

which  is  the  inverse  of  the  transformation  Equation  (A-47). 

We  derive  the  result 


H=  s''P|g‘'s^  +  V 

ijkm 


(A-52) 


by  noting  that  G  ■-  S^g  and  ln^  ^1} ;  where  S  =  •$  i 

H  ^  •  1*1*1 

For  Equation  (A-51)  to  reduce  to  Equation  (A-44)  the  condition 


1  zl 

s^Es"'‘p,s^s-'  =  5^P^ 


I, IT. 


(A-53) 


must  be  fulfilled. 

If  the  kinetic  energy  of  the  rolaling-vibrating  molecule  was  expressed  in  terms  of  the 
Eularian  angles,  0,0,  X  normal  coordinates  Q|.,  together  with  the  conjugate  momenta 

Pe’Pd’  Pz  Pji '  hamiltonian  operator  could  immediately  assume  the  form  of 

Equation  (A-44).  But  the  kinetic  energy  is  expressed  in  lemas  cf 

X  “  P  x)»  (^y  “  P  y) » (^z  ~  Pz)  and  p  .  We  therefore  need  to  apply  the  condition 
Equation  (A-53)  to  the  transformation  of  the  momenta  P  q>  P  P  y  and  P  j.  lo  the  momenta 
found  in  Equation  (A-36). 

To  achieve  this  transformation,  we  first  note  that  the  total  angular  velocity  vector  can 
be  written  as  (see  Figure  A-2) 
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I 


cb  =  9-f  04-  Z  =  Ci)^i  +  0)yj  +  (02k  =  0)j^i'+ COyf  +  OD^k'  (A-54) 
where  the  components  of  3)  are 


0)jj  =  i  &  =  i-  0-fi-  o+i  -  X 


.  Y  ,eic. 


Figure  A-2  -  Components  of  Total  Angular  Velocity .  From  fl)  (p.  2811 


From  Figure  A-2  we  see 


©x  =  sinx  ■  8  -sinScos  X'  <> 


a>y  =  cos  X  -  9  -r  sinGcos  x  •  6 


{A-56) 


0)z=  cos  0  p  -r  X 


with  the  inverse  transformation 


6  =  sinx  ■  (h^-tcos  X’ 


P=:-CSC  0COSX-  fOx-rCSC  GsioX  '  G), 


(A-57) 


X  =  cot  0COS  X  -  (t>x  -  cot  Bsinx '  Oy-r 


Recall  that 


p  =  iL  =  iL  i  ^  £L  ^  =  J2_n  -  -i^D  • 

39^“*  ■  3o^“x  ’  ^“x^8  '  ao^Pp  3e>,Px  (A-5S) 


From  Equation  (A-57)  we  can  find  ^  ,  etc.  resulting  in 


=  sinxpg  -  CSC  Geos  xp©  cot  Geos  xp^ 


^  =  cos  xPq  -i-ese  Gsin  XP©"  6sinxpj^ 


{A-59) 


From  the  definitions  of  px»  Pv»  Pz  i*'  Equation  (A-32)  we  ha^'c 
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P'x  =  (Px-Px)  =sinxpQ-csc  e  cos  xp^  + cote  cos  xp^  -X^kPk 


P'y  =  (Py-py)  =cos  xPe  +  csc  Gsinxp^-cot  Qsinxp^  “S^kPk 
P“z=  (Pz  -  Pr:)  =  P  X  -  kP  k 

Pk=pk 

We  have  now  found  the  coefficients  Smi  of  the  transformation  Equation  (A-51).  The  inverse 
equations  are 


Pq  =  sinxP'  X  +  cos  xP'y  +  x3  k  +  cos  k)  P  k 
p^  =  -sin  0COS  xP'x+  sin0sinxP’y+ cos  0P'z 

+  ^(-  sin  0COS  x3  +  sin0sinx5?  +  cos  0K  j^)Pj^ 

Px  =  P’z+Z^kPk  (A-61) 

Pk  =  Pk 

giving  the  coefficients  s'*^  of  the  inverse  transformation  Equation  (A-47).  Rewriting  the 
condition  (A-53)  as 


I 

ikm 


s"^'p.s 


km 


s“’  =  0 


and  noting  that  pj  is  a  differential  operator  allows  us  to  write 


(A-62) 
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to(p  i®"') + X  '“"C  Pi=  ta)  =  X(p  i®'') 

ikm  ikm  i 


(A-63) 


It  has  been  shown  that  the  above  condition  is  satisfied  for  the  case  in  hand,  allowing  us  to  write 
the  quantum  mechanical  hamiltonian  in  the  form  of  Equation  (A-50) 


1  zl  1  1  zl  1 

+  +V  (A-64) 

gg'  ^ 

As  an  aside,  we  stop  here  for  a  moment  and  quickly  derive  the  hamiltonian  used  as  a 

starting  point  in  the  rigid  rotor-harmonic  oscillator  approximation  presented  in  Chapter  2,  This 

approximation  begins  with  the  hamiltonian  (A-64),  which  is  completely  general  and  contains 

no  approximations.  The  first  approximation  we  make  is  to  neglect  the  dependence  of  |4  and 

l-l  on  the  no"  lal  coordinates  [see  Equation(A-35)I.  This  means  they  are  now  constant,  and  are 

OO 

not  affected  by  any  operators.  Thus  Equation  (A-64)  reduces  to 


H=  iXPjg<Pg  -Pg)(Pg  -P5.)  +  iXp(  + V  (A-64a) 

gg'  k 


Next  we  perform  a  unitary  transformation  to  ensure  the  axes  moving  with  the  molecule 
coincide  with  the  principle  inertial  axes  of  the  molecule  at  equilibrium,  thereby  causing  the  off- 


diagonal  terms  of  the  inertia  tensor  to  vanish.  We  also  neglect  the  terms  of  the  form 


X3,5ik 


since,  by  definition,  they  depend  on  the  squares  of  the  small  vibrational  displacements  (see 
Equation  (A-26)].  Then  the  terms  I'xy  /  ^'yz  and  are  zero,  andrxx=  I  xx'  ^'yy~  ^  yy 
I'zz—  ^zz  tsee  Equation.  (A-32)].  The  terms  kf  gg  now  become  reciprocals  of  the  principle  rigid 
moments  of  inertia  and  the  hamiltonian  is  now  as  shown  in  Equation  (1). 
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A-3  Perturbation  Treatment^o 


A  more  convenient  form  of  Equation  (A-64)  is 


1  11  1-11 


(A-65) 


where 


4 

|I^  |I 


■  P  .P 

gg  ^g 


i  zl] 

.P  .  P 
^g  gg 


The  hamiltonian  in  Equation  (A-65)  can  be  divided  into  three  terms:  Hq  a  purely  vibrational 

TjO 

term  consisting  of  the  last  three  terms  of  Equation  (A-65);  which  consists  of  those  parts  of 


the  first  term  which  are  diagonal  in  the  vibrational  quantum  numbers;  and  the  perturbing 
hamiltonian  H'  which  includes  the  second  term  and  those  parts  of  the  first  term  which  are  not 

included  in  . 

We  now  make  the  following  approximations:  the  wave  function  V  VR  can  be  written  as 
the  product  of  the  two  independent  functions  ¥  y  and  ¥  r  where  ¥  y  are  the  orthonormal 

solutions  to  the  harmonic  oscillator  hamiltonian,  and  ¥  r  depends  only  on  the  rotational 
coordinates.  In  this  manner,  Hq  and  are  diagonal  in  V,  but  not  necessarily  in  R.  If  we 
expand  ¥  r  for  an  asymmetric  rotor  in  terms  of  the  symmetric  rotor  basis  functions 
(see  Appendix  B ),  we  see  ¥r  is  not  diagonal  in  K. 

Before  we  apply  a  perturbation  treatment  to  H  we  must  take  care  of  the  off-diagonal 
elements  of  H'  in  the  vibrational  matrix.  We  do  this  by  performing  a  unitary  transformation 


l^The  theory  of  Section  A-3  is  taken  from  the  following  References:  (3]  (pp.33-37);  [6]  [pp.  263- 
266);  (251  (pp.  277-281). 


96 


which  will  remove  to  first  order  in  X  the  off-diagonal  matrix  elements  of  H',  while  preserving 
to  first  order  in  X  the  matrix  elements  of  H'  that  lie  inside  the  diagonal  blocks  (see  Figure  A-3). 
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Figure  A-3  -  Matrix  Representation  Of  H’.  From  [3]  [p.  34). 


2 

Thus  the  elements  of  H'  which  are  outside  the  diagonal  blocks  are  of  the  order  X  and 
contribute  to  the  energy  only  to  fourth  order  in  X  •  These  elements  can  be  neglected  in  second 
order  perturbation  theory.  This  unitary  transformation  is  known  as  the  Van  Vleck 
transformation. 

A  unitary  operator  can  be  exprcsseo.  as 

U  =  exp(iXS)  (A-66) 

where  S  is  hermitian.  We  transform  H  as  follows 

G  =  U"'‘PiU  =  [l  -  i^s  -  ^  -h  ^  -h  ■  ••  Jho  +  -h  H’j^I  +  iXs-^ - ] 

=  G Q  +  X,G ^-hXG^  +  ^Gj-i--  -  ■  ( A-67) 
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Equating  like  powers  of  X  gives 


Go=Ho+H° 

G,  =  H'+  ((H„  +  H“)S  -  S(H„  +  H°)] 

Gj  =  i(H'S  -  SH')  +  S(H„  +  H°)S  -  +  S^(H„  +  H°)]  (A-68) 

Now  let  the  rotational  sub-block  of  interest  be  labeled  by  |V ;JKM)  and  indexed  by  j,  k, 
1  etc. ,  and  the  levels  of  the  other  rotational  sub-  blocks  by  a  P  T  and  so  forth.  The  Van  Vleck 

transformation  must  satisfy  two  requirements.  Requirement  one  is  that  (|G  Jk)  =(j|H'|k)  or 

0  =  (^((h,,  +  H°js  -  S(H„  +  h;)] k>  =  (E|  -  E  ^)<j|S|k)  (A-69) 

This  requirement  is  satisfied  by 


(jSik>=0 


Requirement  two  is  that  (j^G  j|a)  vanishes,  or 


(A-70) 


0  =  <jH';a>  +  i<j;[(H<,  +  H°  )S  -  S(H„  +  ]a>  =  <jHla)  +  i(E  j  -  E„)  (j|S|a> 

(A-71) 


which  implies 


<|S!a)  = 


KjiH’a) 


E  -E 
J 


a 


(A-72) 


Thus  we  have  determined  the  nature  of  the  hermitian  operator  S  that  generates  the  unitary 
transformation.  Now  by  applying  Equation  (A-70)  and  Equation  (A-72)  to  Equation  (A-67)  we 
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can  determine  the  matrix  elements  of  G  v/ithin  the  same  (2  j  x  1)  x  (2  j  x  1)  block  jV  ;JKM) 


to  second  order 


<jG<,|k>=<j!H„  +  H°|k)=E.8 


<jGJk>  =  <j'H'|k> 


(A-73) 


J  Z 

Now  the  off-diagonal  sub-blocks  contain  terms  only  of  order  A.  and  higher  and  can  be  neglected 
in  second  order  perturbation  theory.  If  more  accuracy  is  desired,  G  can  always  be  subjected  to 
another  Van  Vleck  transformation. 

If  we  now  neglect  the  off-diagonal  elements  of  H  in  the  matrix  V,  the  remaining 
diagonal  terms  each  factor  into  smaller  matrices,  one  for  each  vibrational  state,  whose  elements 
(R,  ViH|R V)  are  labeled  by  the  rotational  quantum  numbers  only. 

Using  standard  perturbation  theory,  we  can  write 


H  =  Ho+Hj^+  H' 


(A-74) 


where  we  can  expand  the  perturbing  hamiltonian  H'  to  second  order  in  a  pov/er  senes  in  X 
( X  small),  that  is 


H’=  rH'2 


(A-75) 


The  zeroeth  and  first  order  terms  arc  given  by 
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(rv|h„  +  H°  +  UijR'V)  =  (R\lH„iR  v>  +  <rv|h°  +  xhIr’ v>  = 


Ev+  iZ<HgsXRMPsP5-  R'V)  -  A.2(hg)<RV|Pg|R'V) 

gg'  8 


(A-76) 


Ev  is  the  vibrational  energy,  diagonal  in  R  and  V,  whose  diagonal  elements  are  Ey  The  last  tv/o 
terms  come  from  the  fact  that  P^gg.  and  hg  are  functions  only  of  the  vibrational  coordinates  (see 

Equation  (A-33)  and  Equation  (A-65)l  plus  the  approximation  V  VR  ~  ^  R  • 

We  can  further  simplify  Equation  (A-76)  by  noting  that,  since  hg  is  pure  imaginary 
(each  operator  pg  carries  one  term  of  -ih  and  all  other  terms  are  real),  and  because  H  and  Pg  are 
hermitian,  the  matrix  ^  is  hermitian.  Also  for  a  non-degenerate  asymmetric  rotor, 

the  V  V  arc  real.  Hence  the  diagonal  term 


<V|hg|V>  =  <v|h'  V)  =  <V|h'g|V>  =  - <V|hj]V>  =  0  (A-77) 

The  second  order  term  is  given  by 

-12^.  _'i2  V  (VR!KlV“R")(VR'iHiVR> 

"2  •"  L  E~rE -  (A-78) 

R  -v"  V  v" 

V"*V 

Here  we  have  assumed  the  rotational  spacing  is  small  compared  to  the  vibrational  levels  so 
that  t  t  can  be  replaced  by  b  ^  -  E  y„ .  Equation  (A-78),  when  multiplied  out  will 
give  us  the  quartic,  cubic  and  quadra  tic  terms  of  ^x,  Py  and  P^.  The  coefficients  of  the  terms 

^  X  ,  A  y  and  A  2  are  of  the  form 


V 


(VI,ggA'-)(V 


ihg|V)^(V,hg.V)(V-HggjV) 


Ev-E, 


V5:V 


(A-79) 
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They  can  be  shown  to  vanish  by  noting  that,  since  hg  is  pure  imaginary, 

(VjhgjV ')  =  -  (VjhgjV") ;  and  since  M-gg  is  a  real  operaU>r,  =  (V  "|M-ggpO  • 

The  rest  of  the  cubic  terms  are  of  the  type  Px  P  y  ~  PyF’x  ,  etc.  These  can  likewise  be 
eliminated  by  use  of  the  commutation  rules  jj  ~  ^  ijk  ^  k  3*  the  cost  of  introducing  extra 
quadratic  temw. 

Thus,  a?  '  result  of  the  perturbation  treatment,  the  rotational  matrix 
H  =  Hg  +  "b  K',  which  corresponds  to  a  single  vibrational  state  can  be  written  as  a 
polynomial  in  P  x ,  P  y  and  P  2 

H  =  E  V  +  Sa  T  S  "  Pj  P”  (A-80) 

83  ®  8  4  ....  88  J)  S  8  J  j 
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88  ;j 


where  and  are  functions  of  the  vibrational  state  only.  Note  also  that 

oc  00 


X  ....=  I 

gs }) 

VstV 


Ev-Ev- 


{A-81) 


A-4  The  Reduced  HamiltoniarJ’ 


The  form  of  the  hamiltonian  in  Equation  (A*80)  has  several  drawbacks.  First,  if  further 

accuracy  is  desired,  higi'.cr  order  perturbation  must  be  considered,  v\  hich  could  become  extremely 

cumbersome.  Also,  not  all  the  coefficients  and  X  are  determinable.  If  all  the 

8S  gg  u 

constants  in  Equation  ( A-SO)  were  known,  it  would  be  a  straightforward  process,  in  principle,  to 
determine  the  energy  levels  from  them.  However,  the  situation  that  we  are  faced  with 
experimentally  is  the  reverse  of  this;  we  wish  to  determine,  from  observed  energy  levels,  the 

l^Tho  theory  of  Section  A-4  is  taken  from  the  following  References:  [4]  Ipp.  .*;-80];  (18J;  [17)  (pp. 
1935-1947). ' 
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values  of  the  rotational  and  centrifugal  constants.  The  problem  in  carrying  out  this  reverse 
calculation  is  that  not  all  constants  or  combination  of  constants  contribute  to  the  energy  levels 
and  are  therefore  not  determinable  from  them.  For  example,  neglecting  the  quartic  distortion 
terms  in  Equation  (A-80)  allows  us  to  write  the  rotational  hamiltonian  as 

^rot  =  2-^%gJgJg-  (A-82) 

gg' 

where  CC  .  are  related  to  the  values  in  Equation  (A-35)  ,  for  the  particular  vibrational  state 

OO 

and  the  symbol  J  now  represents  the  total  angular  momentum.  We  can  reduce  the  matrix  CX  .  to 

OO 

diagonal  form  b\'  a  simple  rotation  of  axes  (to  a  principle  axes  system)  giving  us  the  reduced  form 
of  the  hamiltonian 


h„.=xi^,+  y;^+zj; 


(A-83) 


where  X,  Y,  Z  are  the  principle  values  of  2  ^gg'  •  obvious  that  the  rotational  energies 

depend  only  on  the  principle  values  of  2  ®^gg'  ^nd  all  the  other  coefficients  2  ^gg* 

determinable  from  those  energies.  Therefore,  to  obtain  increased  accuracy  and  to  determine  all 
the  coefficients  in  the  hamiltonian,  it  is  necessary  to  write  the  hamiltonian  in  another  form. 

The  hamiltonian  in  Equation  (A-80)  contains  such  terms  as  x  J  zJ  xJ  y  with  individual 
terms  appearing  in  any  order.  By  using  the  commutation  relations  P  j]  ”  ^ijk  we  can 

rewrite  such  these  terms  in  the  form  2  v  x)  yi  7.  )  zJ  y)  x  J  -i-  terms  to  lower  degree  in  J.  For 

example 


J  xJzJ  xJ  y  ~  2  yJz  J  zJ  yJ  x)  2^  ^  ^  ~  2(^  JzJ  y)  (A-84) 


The  terms  of  low'er  degree  can  then  be  changed  in  the  same  way,  and  the  process  can  be  repeated 
rcsulHng  in  a  sum  of  terms  of  the  following  form 
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(A-85) 


pqr 


where  the  coefficients  hpqr  may  or  may  not  be  complex. 

To  simplify  further,  we  note  that  the  rotational  hamiltonian,  of  a  given  vibrational 
state,  is  invariant  to  both  operations  of  hermirian  conjugation  and  time  reversal.  That  is 


H,o>=H!L.  +  TH„,r'  =  (TH„,T-‘) 


(A-86) 


where  the  symbol  ^  represents  the  adjoint  of  the  operator.  In  addition,  the  operators  J  and  the 
coefficients  hpqr  have  the  follow’ing  properties 

TJT"^  =-]  =  h*pqr  ^  ThpqrT"^  =  h'pqr  <A-87) 

Equation  (A-86)  and  Equation  (A-87)  together  give  us 


pqr 


=  K-  l)'’^‘'"'h*pqr(jxjvj;- 

pqr 


(A-88) 


pqr 


The  last  equation  of  (A-88)  shows  that  since  terms  of  odd  degree  in  p+q-rr  change  sign  under 
operations  of  time  reversal  and  hermitian  conjugation,  they  do  not  satisfy  the  invariance 
relation  (.'\-S6).  The  first  two  equations  of  (A-S8)  show  that  the  remaining  terms  of  even  degree 
satisfy  the  invariance  rcLlion  only  if  the  hpqr  are  real.  We  now  have  a  simplified  rotational 
hamiltonian 
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(A-89) 


p+q-fr 
even 

where  hpqr  are  real.  Note  that  the  coefficients  of  the  terms  of  degree  greater  than  two  are  the 

centrifugal  distortion  coefficients  for  a  particular  vibrational  state  of  the  molecule.  One  further 
simplification  can  be  made  by  noting  that  for  molecules  of  C  (orthorhombic)  symmetry  the 

only  non-zero  terms  in  Equation  (A-89)  are  those  with  p,  q,  and  r  all  even,  giving  us 

^  rot  ~  X  h  pqr  (j  xJ  yJ  z  J  J  yJ  x)  'hpqr  real  (A-90) 

P'q'f 
even 

If  further  accuracy  is  desired,  wo  can  just  extend  the  power  series  to  any  degree  in  J  that 
we  want.  For  the  purposes  of  this  experiment,  expansion  of  the  power  senes  (A-90)  to  the  sixth 
power  in  J  is  needed.  With  a  slight  change  in  notation  we  have 


a?  1  7  7  7  V  \ 

x^ry-^j;+j;jyrx) 


(A-91) 


where  the  coefficients  are  real  and  ^  jjp  ”  terms  of  J  ,  J  z  and  J  ±  we  have 


Hrot  = 


020- 


220 


040- 


6C0 


4^ 


^202 '*'*^022 *^420  ^  "^^222^ 


where  represents  the  anti<ommutator  AB  -f  BA  and  the  subscripts  on  the 

coefficients  now  refer  to  the  powers  of  J  .  J  z  and  J  r  respectively.  To  go  from  the  form  of  the 
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Hrot  Equation  (A-91)  to  Hrot  Equation  (A-92)  requires  repeated  use  of  the  conunutation 
relations  jj  j' J  jj  ~  ^  ijk  J  k  •  For  each  constant,  terms  of  higher  degree  are  negligible  therefore 

neglected.  With  this  in  mind  the  relations  between  the  constants  in  Equation  {A-91)  and 
Equation  (A-92)  are  given  in  Table  A-1. 
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■ri^r ^  -w  ‘ 


TABLE  A-1 

Relationship  among  constants  in  Equation  (A-91)  and  EquaHon  (A-92).  From  [4]  [p.  21]. 


®  200  ~  2(^  X  ®y) 


sC^xxz+^yyz  ^xyz)"*^; 


^002  4(®x  By)  ^006  “  64(^»«  “^Wy)“‘^(^xxy -<^yyx) 

^20=^^-^00-^6T^ 

"^400  ~  ^^xy) 

^220  ~  C^xz  ■*■  ”^>2)  ~  ^”^400 

T  =  T  -  T  -  T 
040  ®  220  ^  400 

^202  “  4(^>«  “^yy) 

^600  ”  16  XXX  +  ^yyy)  +  sC^xxy  ■*'^yyx) 

^420  “  4(^  xxz  +  ^yyz)  +  i^xyz  “  ^^eOO 
^240  (^22^  "^^2^  )~  ^240  “^^600 

^060  ^  ^222  -  ^240  “^420  ~  ^600 

^402  “  64  (^XXX  ”  ^  yyy)  +  xxy  ~  ^yyx) 

^222  =i-(^  xxz 20,02 

^042  ~  ^  ~  )  “  ^222  “  ^402 

^204  “  ■^(‘^xxx  ^yyy)  ~  ^(^xxy  +^yyx) 


We  still  have  the  problem  of  which  constants  or  combination  of  constants  contribute  to 
the  energy  levels  and  can  therefore  be  determined  from  them.  We  can  solve  this  problem  by 
noting  that  the  eigenvalues  of  the  rotational  hamiltonian  are  unaltered  when  the  hamiltonian 
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is  transformed  by  an  arbitrary  unitary  operator  [the  transformation  from  Equation  (A-82)  to 
Equation  (A-83)  is  an  example  of  such  a  unitary  transformation].  If  we  construct  the  unitary 
operator  as  a  power  series  in  the  total  angular  momentum  J  and  other  parameters,  then  the 
transformed  hamiltonian  will  be  a  power  series  similar  to  the  original  hamiltonian,  but  with 
coefficients  depending  on  the  parameters  in  the  unitary  transformation.  Because  these 
parameters  are  arbitrary,  we  can  choose  the  unitary  operator  in  such  a  way  as  to  eliminate  as 
many  terms  as  possible  from  the  transformed  hamiltonian.  The  coefficients  of  the  remaining 
terms  are  the  maximum  number  of  terms  which  can  be  determined  from  the  experimental  data. 
The  situation  is  not  quite  as  simple  as  described  above,  but  we  can  still  make  a  meaningful 
transformation  by  considering  orders  of  magnitudes  of  various  coefficients  involved,  and 
discarding  terms  which  may  be  negligible.  Details  of  this  process  can  be  found  in  Reference  [4]. 

A  hamiltonian  subjected  thus  to  a  unitary  operation  is  called  a  reduced  hamiltonian. 

We  wish  to  construct  a  reduced  hamiltonian  Hf  with  the  same  eigenvalues  as  . 

The  reduced  hamiltonian  is  given  by 

H'  =UH  U“’  (A-93) 

rot  rot 

Recall  that  the  most  convenient  form  for  a  unitary  operator  is 

U=exp(iS)  (A-94) 

/ 

where  S  is  hermitian.  To  make  H  invariant  under  time  reversal,  we  must  make  U  invariant, 
which  requires  that  S  change  sign  under  time  reversal.  Therefore,  if  we  express  S  as  a  power 
series  similar  in  form  to  Equation  (A-85),  we  see  it  has  real  coefficients  and  contains  powers  of 
p+q+r=odd  only. 


(A-95) 


p+q+r 
odd 
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The  exact  stracturc  of  S,  and  therefore  of  U,depends  on  the  symmetry  of  the  molecule  in  question. 
IDetails  of  the  arguments  and  exact  forms  of  U  and  S  can  be  found  in  Reference  [4]. 

As  a  result  of  these  arguments  it  has  been  shown  that  for  molecules  of  orthorhombic 
symmetry,  the  hamiltonian  of  Equation  (A-91)  can  be  transformed  to  an  analogous  form  by 

the  following  changes  in  coefficients 


/ 

nn  —  nr 

^  XX  “  ^  XX 

T  xy  =  Txy  -2(B  X- 

^  XXX  ~  ^xxx 

^  )0y  “  ^xxy  “  X  “  ®  y)S3ii  “  4(Txx  -  ^xy)s^ll  ~  4(B  x  “  B  (A-96) 

^'xx2='I'xxz  +  2(Bx  -B,)S3j^  +4(T^  -4(Bx -6^)52^^ 

^  xyE~  *^xyz"'^[(^y“  ^z)^311  '^(®z“^x)^131  ■*'(®x~  ^  7)^113] 
where  S  S  i  ,  and  S  ^^3  are  coefficients  from  Equation  (A-95) 

H  of  Equation  (A-92)  can  also  be  transformed  to  an  analogous  form  H  by  making 
the  coefficient  transformations  indicated  in  Table  A-2.  Note  the  coefficients  S  , 

and  S  have  been  replaced  by  the  parameters  p ,  X,  and  V  which  are  also  defined 

in  Table  A-2. 
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TABLE  A-2 

Coefficients  of  the  Transformed  Hamiltonian.  From  [4]  [p.  33]. 


®  200 


T 

^  400 


T’  = 

202 


=  B  +  4o 

200  ^ 

^  020  “  ®  020  ~ 

®  002  ®  002 

Q. 

CN 

1 

II 

220  ~  '^220 

040  ~  040  P 

=  T 

202 

022  ""^^022 

Cl 

+ 

II 

Tf 

o 

o 

H 

CD' 


600  = 


<D  -X 

600 


0)'  = 
420 

O'  = 
240 

CD'  = 
^060 


^420  + 

CD -5^  =  104 

^060  “ 


CD'  = 
402 


<D  —  V 

402 


O'  = 

^  222 


O  -  2aX  +  lOv  + 

222 


('-2T^^^ +4oT^_ -8T  -160^0 

V  220  202  004 

B 


O'  = 

^  f\A*> 


002 

[-4T„.„+4aT _ +8T. 


042  ^042 


4>_  -  20H  -  9V+  0«  -022  “'004  '°^ 


B 


002 


O'  = 

204 


O  +  - 

204  2 


O'  = 

^024 


^002 


O'  = 

^006 


^006 


^  ^002^111 


^‘  =  ®  TO  (2^,3  -  ®  31,  -  ^3, )  +  S  -  2B^  s2„ 


^  -  B  /'s  +  2T^„,  s 


002  V  311 


202  111 


''  4^  002(^311  ^  131 )  ■^^^004  ^111  ®  002  ^111 


B 


cy  = 


020 


B 


002 
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The  reduction  of  the  rotational  hamiltonian  is  completed  by  making  appropriate 
reductions  of  certain  coefficients  which  will  fix  the  parameters  p ,  A,,  M-  and  V.  Because 

SF2O2  approaches  the  symmetric  prolate  limit,  we  need  to  consider  a  reduction  that  stays  finite 
in  this  limit.  The  choice  we  make  is  called  the  asymmetric  top  reduction  or  the  S-reduction. 
This  reduction  is  obtained  by  the  conditions 


022 


022 


042 


024 


=  0 


These  conditions  are  solved  to  give  the  values  of  p  ,X,  .d-  and  V  given  in  Table  A-3. 
The  form  of  the  S-reduced  hamiltonian  is 


(A-97) 


H -  D  j(l^)  V,  -  D 

+  d  jJ^Oi  +  J^-)  +  d  j(jt  +  jl)  +  H,( 

+  H^j^  +  h,(rt\ji  +  ji)  +h2j^(jt  +  jl)  +h3(jt  +  j-)  (A-98) 

The  relationship  of  the  coefficients  in  Equation  (A-98)  to  coefficients  from  previous  expressions 
is  given  in  Table  A-3. 
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TABLE  A-3 


Constants  of  the  S-Reduced  Haniiltonian.  From  [4]  [p.  35]. 


-Dj=To04-i°''T„2  -‘=JK=T22„  +3a-'T^  -Dk=T(^ 


-a"^T 

2^  i  022 


:T 

202 

<12  = 

”^004  ■*' 

-a'^T 

4  022 

H  =cE) 

J  600 

-X 

H,k 

=  0 

420 

+  6A,  ~ 

■3]i 

Hk, 

~  ^240 

-5A.+  10|i 

~  ^060 

-7\i 

N  = 

:  0 

402 

■  V 

‘’2  = 

^204 

it 

2 

^3“^006 

+  V 

p= 

-a“^T 

4^^  ^022 

p  = 

(Za^d- 

!)■'{ 

*^042 

024 

-5aT, 

v=  — a"V+ -cj'^0  +  — 

16  ^  8  024  B 


T  T 
004  022 


X  5a  +  2  +  ^222  "^  (  2^  ^220  "'''^202  ^  ^022  ^”^004) 
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I 

I 

I 

I 

I 

I 

I 

f 

I 

i 

I 

I 

I 

I 

I 

I 

I 

I 
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Appendix  B 

Wave  Equation  Solution  For  The  Symmetric 

Topi2 

The  motion  of  a  symmetric  top  is  usually  described  by  the  Euler  angles  (see  Figure  B-1). 

z 


Figure  B-1  -The  Euler  angles.  From  (25]  (p.  78  ]. 


Now  let  r  be  an  arbitrary  vector  having  Cartesian  compone.nis  ~  (^X'^Y'^z) 
space-fixed  (lab)  frame  and  the  Cartesian  components  fg  =  (Tx/l" in  the  body-fixed 

^^The  theory  of  Appendix  B  is  taken  from  the  following  references;  (2)  (pp.  60-61);  [25]  (pp.77-83l 
(l][pp.  285-286). 


112 


I 

I 

I 

i 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

s 

I 


(molecular)  frame.  Since  we  are  describing  the  rotation  of  the  sym/.-ietric  top  by  the  Euler  angles, 
the  two  representations  of  the  vector  r  are  related  by  the  unitary  transformation  <1> 


^F-2^Fg(G,o,x)’^g 

O 


^8~  Sr^gF(0,b,x)^F 

F 


where  O  is  called  the  direction  cosine  matrix.  Naturally  we  express  the 
bv  the  product  of  the  three  Euler  angle  rotations  (see  Figure  B-2) 


^z(X)^N(e)^Z(0) 


(B-1) 

opeiL-tor  ^ 


(B-2) 
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Figure  B-.?  -  Components  of  Euler  Angle  Rotations.  From  {25]  Ip.  79). 


Since  we  can  reprojent  the  Eule.  '.ngles  as  the  product  of  three  rotations,  we  can  find  the 
relationship  between  Fpand  fg  to  be  as  depicted  in  Table  B-1.  Because  the  are  the  real 


elements  of  a  unitary  transformation. 


That  means  the  transformation  is 


orthonormal,  i.e. 
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Figure  B-2  -  Components  of  Euler  Angle  Rotations.  From  [25J  Ip.  791. 


Since  we  can  reprojent  the  Eule.  "ingles  as  the  product  of  three  rotations,  we  can  find  the 
relationship  between  Tpand  fg  to  be  as  depicted  in  Table  B-1.  Because  the  are  the  real 


elements  of  a  unitary  transformation. 


That  means  the  transformation  is 


orthonormal,  i.e. 


114 


(B-3) 


g 


Fg 


=  5 


FF- 


X^Fg^Fg-  ^gg 

and  Table  B-1  can  be  read  either  across  or  down.  The  entries  in  Table  B-1  are  the  cosines  of  the 
angles  between  the  various  pairs  of  axes. 


Table  B-1 

Direction  Cosines  Relating  Rotating  to  Nonrotating  Axes  as  Functions  Of  Eularian  Angles. 
_ Fromflj  [p.2861. _ _ 


X 

Y 

Z 

X 

cos  0  cos  (?  cos  X 

-sindsinx 

cos  0 Sind  cos  x 

+  cos  dsinx 

-  sin  0  cos  X 

y 

1 

1 

-cos  8  cos  0  sin  x!  -  cos  0sin<|)sinx 

-Sind cos  X  +  dcos  x 

i 

sinOsinx 

I 

sin  0  COSO  ;  sin0sin({) 

cos  0 

By  examining  Figure  B-1  we  deduce  the  following  relations 


A 


A 

=  n 


z 


ng=-n^sin4)  +  nYCOS(i) 


(B-4) 


fi.  =  n  sin0cos  0  +  n  sin0sin(|>  +  n  cos  0 

XX  • 

By  noting  that  the  projection  of  the  angular  momenU'm  J  on  an  arbitrary  axis  n  is  the  operator 
which  generates  an  infinitesimal  rotation  about  ii,  that  is 
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(B-3) 


=  8 


FF- 


S^Fg^Fg-  ^gg- 

and  Table  B-1  can  be  read  either  across  or  down.  The  entries  in  Table  B-1  are  the  cosines  of  the 
angles  between  the  various  pairs  of  axes. 


Table  B-1 

Direction  Cosines  Relating  Rotating  to  Nonrotating  Axes  as  Functions  Of  Eularian  Angles. 


From  (iHp.  286). 

X 

Y 

z 

X 

cos  0  COS  <)  cos  X 

-sin(i)sinx 

cos  Gsintj)  cos  x 

+  cos  (|)sinx 

-sinGcos  X 

y 

-  cos  0COS  <?sin  % 

-sin(})cos  y 

-  cos  0sin(})sinx 

+  cos  (j)COS  X 

i 

sinGsinx 

z 

sinGcos  0 

sin  9  sin  0 

cos  0 

By  examining  Figure  B-1  we  deduce  the  following  relations 


n  =-n  sin4>  +  n  cos 

0  X  Y 


0 


{B-4) 


fi  =  n  sinGcos  0  +  n  sin  9 sin 6  +  n  cos  0 

^  A  Y  ^ 

By  noting  that  the  projection  of  the  angular  momenv.»m  J  on  an  arbitrary  axis  n  is  the  operator 
which  generates  an  infinitesimal  rotation  Ctp  about  n,  that  is 
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h 


h 


h 


(B-10) 


A  = 


Sr.^l , 


B  = 


Sit^l 


C  = 


Equation  (B-10)  implies 


A>B>C  (B-11) 

For  a  symmetric  prolate  top  we  identify  the  symmetry  axis  (z-axis)  as  the  a-axis,  and  the  other 
two  axes  as  the  b-axes.  This  gives  us  the  hamiltonian 


H=  BJ^-f  J^(A -B)  (B-12) 

We  use  the  hamiltonian  (B-12)  in  the  Schroedinger  equation.  We  can  calculate  J  and  J  2 
from  Equation  (B-7) 


_J _ ^ 

sin  6  96 


sin  ^6  96^ 


2cose  9“  f  1  \9^ 

sin  ^6  sin  ^0  /  9x^ 


Upon  substitution  into  Schroedinger's  equation  we  get 


(B-IS) 


-sin  ^0 


^(sin6i)+^ 


2  2 
sin  0 


cos^0  A 

-  2„  B 
sin  0 


2cos  0  9 

.in  20  ^r,dO 


0  {B-14) 


We  assume  a  solution  of  the  form 


(B-15) 


where  M  anc'  K  must  be  integers  to  preserve  the  periocity  of  the  solution.  0^^^  must  now  satisfy 


1 17 


the  equation 


V  sine)  d0 


By  letting 


) 

f  cos  h 

aI 

^2cos  0'\ 

E 

sin  e 

\  sin  ^0 

^BJ 

F  +( 

V.  sin  ^0  J 

KM  + 

{B-16) 


we  get 


x  =  ^(1-  cos  0) 


(B-17) 


e,„>=x^  (1-x)-  F 


(x) 


(B-18) 


where  is  a  polynomial  in  x.We  find  the  equation  for  to  be 


x(l  -x)f"^^)+  (a-px)F'(y)+YF(^j  =  0  (B-19) 


where 


a  =  iK  -  M!  -5-  i 


P  =  !K  +  Mi  +  iK  -  iMI  +  2 


(B-20) 


Y=  -  ^  +  K‘  -  (|1K  +  M;  +  -  Mf)(|iK  +  Mi  +  ^iK  -  Mi  + 1) 

«o 

Equation  (B-19)  can  be  solved  by  normal  power  scries  methods.  By  lotting  jj)  =  n 

n=0 

get  the  recurrence  relation 

r  nCn+D-j-pn-yl 

(8-2!) 

n+1  L  (n+l)(n+o)  J  " 

For  a  physical  solution  the  recurrence  relation  must  terminate  which  means  we  get  for  the  energy 
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{B-22) 


E  =  BJ(J  -f  1)  +  (A  -  B)K^ 

where 

J  =  nma.  "vlK  +  Ml+ilK-Ml  (B-23) 

and  ri  max  is  the  largest  value  of  the  integer  n  for  wh'ch  the  recurrence  relation  (B-21)  does  net 
vanish.  Note  from  Equation  (B-23)  that  J  (the  total  angular  momentum)  must  be  a  positive 
integer  such  that 


J  =  0,1,2,... 

K=0,±  i,±2,...,±  J  (B-24) 

M=04: 1,±  2,...,±J 

Rnally ,  note  that  K  and  M  represent  the  projection  of  the  total  angular  momentum  J  onto  the 
space-fixed  (lab)  axis  and  the  molecule-fixed  Z  axis  (axis  of  symmetry)  respectively. 
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